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Abstract 

The manifestly gauge invariant exact renormalisation group provides a frame- 
work for performing continuum computations in SU(N) Yang-Mills theory, 
without fixing the gauge. We use this formalism to compute the two-loop (3 
function in a manifestly gauge invariant way, and without specifying the details 
of the regularisation scheme. 
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1 Introduction 

In ref. [1], a manifestly gauge invariant exact renormalisation group (ERG) was in- 
troduced for SU (N) Yang-Mills theory that can be straightforwardly renormalised 
to any loop order. As such, this formalism is suitable for general continuum calcu- 
lations, which can be performed without fixing the gauge. The obvious obstruction 
to not fixing the gauge — namely that the inverse of the gauge field two-point vertex 
does not exist — is avoided by the particular choice of quantities we aim to compute 
(just as on the lattice). In this, and previous works [2-16], we compute vertices of 
the Wilsonian effective action; in the future we aim to generalise these methods to 
compute the correlators of gauge invariant operators [17]. 

The construction of a real, gauge invariant cutoff, A, is achieved by embedding the 
physical SU(N) gauge theory in a spontaneously broken SU(N\N) gauge theory [9]. 
To compute the effective action without fixing the gauge, we use the fact that there 
an infinity of possible ERGs that specify its flow as modes are integrated out [18] (the 
continuum equivalent to the infinite number of ways of blocking on a lattice [1, 8]) 
and that out of these there are infinite number that manifestly preserve the gauge 
invariance. We then further specialise to those ERGs (still infinite in number) which 
conveniently allow renormalisation to any loop order [1, 16]. 

The key to doing this is first identifying which of the infinity of dimensionless 
couplings in the regularised theory must be renormalised. For technical reasons, the 
superscalar field which spontaneously breaks the SU(N\N) symmetry is given zero 
mass dimension [13], and thus is associated by the usual dimensional reasoning with 
an infinite number of dimensionless couplings. That these couplings do not require 
renormalisation has been assumed in [1] but will be proven in sec. 3.5 (see also [16]). 
Thus, the only couplings requiring renormalisation are <?(A), which is associated with 
the physical SU(N) theory, and 52(A), which is associated with an unphysical copy 
present due to the SU(N\N) regulating structure. For convenience, it is useful to 
define a quantity 

a:=g 2 2 /g 2 . (1) 

To consider on-shell gluons, one can gauge fix after the computation [13]. 
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Since the manifest preservation of gauge invariance ensures that the gauge fields have 
no wavefunction renormalisation [3], g and 52 are the only quantities which run. 

To facilitate calculations in our manifestly gauge invariant ERG, a powerful set 
of diagrammatic techniques have been developed [1, 16, 19]. In this paper we com- 
plete the description of the associated computational scheme and demonstrate its 
consistency and potential by computing the SU(N) Yang-Mills two-loop (5 function, 
02, without fixing the gauge. 

At the heart of the diagrammatic calculus is an elegant representation of the 
flow equation. Its simplicity arises, counterintuitively, from the immense freedom 
in the precise construction of the formalism. We are able to turn this freedom to 
our advantage by recognising that many of the details of the setup are non-universal 
and, moreover, need never be explicitly defined. Knowing that such details must 
necessarily cancel out in the computation of a universal quantity, we can efficiently 
absorb them into diagrammatic rules. Thus, the diagrammatic flow equation hides 
a terrific amount; whilst these hidden features must be properly understood when 
constructing the formalism [1, 16], they can be essentially forgotten about when it 
comes to performing actual calculations [19, 20]. 

Within our ERG, the flow is controlled by a (generically) non-universal object, 
S, the 'seed action' [1, 10, 13, 14, 16,21]. This respects the same symmetries as the 
Wilsonian effective action, S, and has the same structure. However, whereas our 
aim is to solve the flow for S, S acts as an input. By choosing the two-point, tree 
level seed action vertices equal to their Wilsonian effective action counterparts, we 
can arrange for the ERG kernels, integrated with respect to din A, to be the inverses 
of the classical two-point vertices, in the transverse space [1, 13, 14, 16]; equivalently 
we can say that they are inverses up to remainder terms. These are called 'gauge 
remainders' and exist because the manifest gauge invariance demands that they are 
there. In recognition of both their role and form, we refer to the integrated ERG 
kernels as effective propagators, mindful that they are not propagators in the usual 
sense. 

The diagrammatic procedure for computing is as follows. We start by using 
the flow equations to compute the flow of the two-point vertex corresponding to 
the physical SU (N) gauge field, which we suppose carries momentum p. To obtain 
a solvable equation for we specialise to the appropriate loop order and work 
at 0(p 2 ); this latter step constrains the equations by allowing the renormalisation 
condition for the physical coupling <?(A) to feed in. 

We now recognise that certain diagrams generated by the flow comprise exclu- 
sively Wilsonian effective action vertices joined together with an effective propagator 
struck by —Ad^\ a . These terms are processed by moving the —Ad^\ a from the effec- 
tive propagator to the diagram as a whole, minus correction terms in which — Ac^\| a 
strikes the vertices. The former diagrams are called A-derivative terms; the latter 
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can be processed using the flow equation. Amongst the terms thus generated are 
those which can be simplified by applying the effective propagator relation. Such 
terms cancel non-universal contributions up to gauge remainders which can, them- 
selves, be processed diagrammatically. Iterating the diagrammatic procedure, the 
expression for /3 2 ultimately reduces to the following sets of diagrams: 

1. A-derivative terms; 

2. 'a-terms', consisting of diagrams struck by d/da; 

3. 'C(p 2 )-terms', which comprise an 0(p 2 ) stub i.e. a diagrammatic component 
which is manifestly 0{p 2 ). 

The 0(p 2 )-terms can be manipulated. In the calculation of at any rate, the 
structure attaching to the stub can be directly Taylor expanded to zeroth order 
in p — which can once again be done diagrammatically. The above diagrammatic 
procedure is then iterated. At two loops, as mentioned in [1], this procedure is 
not so straightforward, since naive Taylor expansion can generate infra-red (IR) 
divergences. 

The strategy for dealing with such diagrams is to recognise that, by considering 
sets of terms together, these IR divergences cancel out. Organising the calculation 
in this way is facilitated by diagrammatic 'subtraction techniques', which we de- 
scribe in sec. 3.4. Now the 0(p 2 ) can be processed, and so fa can be reduced to 
just A-derivative and a-terms. As anticipated in [1], agreement of with the stan- 
dard, universal answer is expected only in the limit that a — > 0. In section 4.5 we 
demonstrate that, subject to some very general constraints, the a-terms vanish, in 
this limit. 

It is from the A-derivative terms that the universal coefficient can be extracted, 
which is most easily done by working in D = 4 — 2e. At the one-loop level, life is easy: 
each individual diagram is finite; the leading order contribution to (3\ is finite and 
universal, with all non-universal contributions vanishing as D — > 4. At two-loops, as 
one would expect, individual diagrams can develop divergences as D — ► 4. Though 
these must of course cancel between terms, the surviving finite contributions do not 
obviously combine to give something universal. The trick is to once again employ the 
subtraction techniques. This allows us to isolate non-universal contributions, which 
then cancel amongst themselves. The remaining contributions can be evaluated 
directly, combining to yield the expected answer. 

This paper is organised as follows. In sec. 2 we review the setup and the various 
diagrammatic techniques that we will require. In sec. 3 we describe the technology 
for evaluating the A-derivative terms. Following a statement of the basic idea, the 
principles are illustrated in the context of a computation of f3\. The general con- 
siderations for two-loop integrals are discussed and then the subtraction techniques 
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are explained. Finally, we give the proof that all dimensionless couplings besides g 
and 52 can be prevented from running. In sec. 4 we give the A-derivative and en- 
terals that contribute to 02 i n the D — > 4 limit and extract the universal, numerical 
coefficient. We conclude in sec. 5. 

2 Review 

2.1 Elements of SU(N\N) Gauge Theory 

We regularise SU(N) Yang-Mills theory by instead working with SU(N\N) Yang- 
Mills theory. The gauge field is valued in the Lie superalgebra and thus takes the 
form of a Hermitian supertraceless supermatrix: 



Here, A l ^{x) = A^rf is the physical SU(N) gauge field, rf being the SU(N) gener- 
ators orthonormalised to tr(rf rf) = S ab /2, while A 2 ^(x) = A\ t% is a second unphys- 
ical SU(N) gauge field. When labelling e.g. vertex coefficient functions, we often 
abbreviate A 1,2 to just 1,2. The B fields are fermionic gauge fields which will gain a 
mass of order A from the spontaneous breaking; they play the role of gauge invariant 
Pauli-Villars (PV) fields, furnishing the necessary extra regularisation to supplement 
the covariant higher derivatives. 

To unambiguously define contributions which are finite only by virtue of the 
PV regularisation, a preregulator must be used in D = 4 [9]. This amounts to 
a prescription for discarding otherwise non-vanishing surface terms which can be 
generated by shifting loop momenta; we use dimensional regularisation. 

The theory is subject to the local invariance: 



The first term, in which = d^ — iA^, generates supergauge transformations. Note 
that the coupling, g, has been scaled out of this definition. It is worth doing this: 
since we do not gauge fix, the exact preservation of (2) means that none of the fields 
suffer wavefunction renormalisation, even in the broken phase [13]. 

The second term in (2) divides out the centre of the algebra. This 'no A shift 
symmetry' ensures that nothing depends on .A and that A has no degrees of free- 
dom. We adopt a prescription whereby we can effectively ignore the field A , alto- 
gether, using it to map us into a particular diagrammatic picture [1, 16]. 




<u M = [v ll ,n(x)} + \^{x)i. 



(2) 
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The spontaneous breaking is carried by a superscalar field 




C 

which transforms covariantly: 

6C = -i[C,0\. (3) 

It can be shown that, at the classical level, the spontaneous breaking scale (ef- 
fectively the mass of B) tracks the covariant higher derivative effective cutoff scale 
A, if C is made dimensionless (by using powers of A) and S has the minimum of its 
effective potential at: 

1 



< c > = "-(o -ij- < 4 > 

In this case the classical action So also has a minimum at (4). At the quantum level 
this can be imposed as a constraint on S, which can be satisfied by a suitable choice 
of S [13, 16]. When we shift to the broken phase, D becomes a super-Goldstone 
mode (eaten by B in unitary gauge) whilst the C l are Higgs bosons and can be 
given a running mass of order A [3, 9, 13]. Working in our manifestly gauge invariant 
formalism, B and D gauge transform into each other; in recognition of this, we 
define the composite fields Fm = (-B«,-D), Fjy = (B u ,—D), where M, N are five- 
indices [1, 16]. 2 

The renormalisation conditions for the couplings g and §2 are: 

S [A = A\C = a] = ^ str jd D x (i^)' + • • • , (5) 



S[A = A 2 ,C = a] = ^^Jd D x (F%) +•••, (6) 



where the ellipses stand for higher dimension operators and the ignored vacuum 
energy. 

2.2 The Flow Equation 

The flow equation is most naturally phrased in its diagrammatic form, as shown 
in fig. 1. As in our previous works, we have not drawn the improperly regulated 
diagram in which the kernel bites its own tail, having removed such terms by placing 
suitable constraints on the covariantisation [1,5, 13]. 



2 The summation convention for these indices is that we take each product of components to 
contribute with unit weight. 
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Figure 1: A diagrammatic representation of the flow equation. 



The left-hand side depicts the flow of all independent Wilsonian effective ac- 
tion vertex coefficient functions, which correspond to the set of fields, {/}. Each 
coefficient function has associated with it an implied supertrace structure (and sym- 
metry factor which, as one would want, does not appear in the diagrammatics) . For 
example, 




represents both the coefficient functions S clcl and S cl ' Cl which, respectively, are 
associated with the supertrace structures strC^C 1 and strC 1 strC 1 . 

The objects on the right-hand side of fig. 1 have two different types of component. 
The lobes represent vertices of action functionals, where £ 5 = g 2 S — 2S. The object 
attaching to the various lobes, * is the sum over vertices of the covariantised 
ERG kernels [4, 13] and, like the action vertices, can be decorated by fields belonging 
to {/}. The fields of the action vertex (vertices) to which the vertices of the kernels 
attach act as labels for the ERG kernels though, in certain circumstances, the par- 
ticular decorations of the kernel are required for unambiguous identification [1, 16]. 
However, in actual calculations, these non-universal details are irrelevant. We loosely 
refer to both individual and summed over vertices of the kernels simply as a kernel. 
Note that we restrict the choice of kernels such that those labelled at one end by 
either A or B and at the other by either C or D do not exist [13]. 

The rule for decorating the complete diagrams on the right-hand side is simple: 
the set of fields, {/}, are distributed in all independent ways between the component 
objects of each diagram. 

Embedded within the diagrammatic rules is a prescription for evaluating the 
group theory factors. Suppose that we wish to focus on the flow of a particular 
vertex coefficient function, which necessarily has a unique supertrace structure. On 
the left-hand side of the flow equation, we can imagine splitting the lobe up into a 
number lobes equal to the number of supertraces (traditionally, these lobes would be 
joined by dotted lines, to indicate that they are part of the same vertex [1, 16]). To 
finish the specification of the vertex coefficient function, the lobes should be explicitly 
decorated by the fields, {/}, where the fields are, in this picture, to be read off each 



S 



lobe in the counterclockwise sense. 

On the right-hand side of the flow equation, things are slightly more complicated. 
Each lobe is, in principle, a multi-supertrace object and the kernel can attach to any 
of these supertraces. The kernel itself is a multi-supertrace object — for more details 
see [1, 16] — but, for our purposes, we need note only that the kernel is, in this more 
explicit diagrammatic picture, a double sided entity. Thus, whilst the dumbbell like 
term of fig. 1 has at least one associated supertrace, the next diagram has at least 
two, on a account of the loop. If a closed circuit formed by a kernel is devoid of fields 
then it contributes a factor of ±iV, depending on the flavours of the fields to which 
the kernel forming the loop attaches. This is most easily appreciated by defining the 
projectors 

and noting that stru± = ±iV. In the counterclockwise sense, a a + can always be 
inserted for free after an A 1 , C 1 or F, whereas a ct_ can always be inserted for free 
after an A 2 , C 2 or F. 

The above prescription for evaluating the group theory factors receives 1 /N cor- 
rections in the A 1 and A 2 sectors. If a kernel attaches to an A 1 or A 2 , it comprises a 
direct attachment and an indirect attachment. In the former case, one supertrace as- 
sociated with some vertex coefficient function is 'broken open' by an end of a kernel: 
the fields on this supertrace and the single supertrace component of the kernel are 
on the same circuit. In the latter case, the kernel does not break anything open and 
so the two sides of the kernel pinch together at the end associated with the indirect 
attachment. This is illustrated in in fig. 2; for more detail, see [1, 16]. 




Figure 2: The 1/N corrections to the group theory factors. 



We can thus consider the diagram on the left-hand side as having been unpack- 
aged, to give the terms on the right-hand side. The dotted lines in the diagrams with 
indirect attachments serve to remind us where the loose end of the kernel attaches 
in the parent diagram. 
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2.3 Ward Identities 



All vertices, whether they belong to either of the actions or to the covariantised ker- 
nels are subject to Ward Identities which, due to the manifest gauge invariance (2,3), 
take a particularly simple form. In this paper, we need only the Ward identities for 
the action vertices, which are illustrated in fig. 3, though note that the Ward iden- 
tities for the vertices of the kernels are very similar [4, 5, 13, 16, 19]. 



x l P v X i Y X Y X . Y X Y 

Figure 3: The Ward identities for action vertices. 



On the left-hand side, we contract a vertex with the momentum of the field which 
carries p. This field — which we will call the active field — can be either A p , A 2 p , Fr 
or Fr. In the first two cases, the open triangle > represents p p whereas, in the 
latter two cases, it represents pr = (p p , 2) [1, 16]. (Given that we often sum over all 
possible fields, we can take the Feynman rule for > in the C-sector to be null.) 

On the right-hand side, we push the contracted momentum forward onto the field 
which directly follows the active field, in the counterclockwise sense, and pull back 
(with a minus sign) onto the field which directly precedes the active field. Since 
our diagrammatics is permutation symmetric, the struck field — which we will call 
the target field — can be either X, Y or any of the un-drawn fields represented by 
the ellipsis. Any field(s) besides the active field and the target field will be called 
spectators. 

The momentum routing follows in obvious manner: for example, in the first 
diagram on the right-hand side, momentum q+p now flows into the vertex. In the 
case that the active field is fermionic, the field pushed forward / pulled back onto is 
transformed into its opposite statistic partner. There are some signs associated with 
this in the C and -D-sectors, which we will not require here [1, 16]. 

The half arrow which terminates the pushed forward / pulled back active field 
is of no significance and can go on either side of the active field line. It is necessary 
to keep the active field line — even though the active field is no longer part of the 
vertex — in order that we can unambiguously deduce flavour changes and momentum 
routing, without reference to the parent diagram. 
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2.4 Taylor Expansion of Vertices 



For the formalism to be properly defined, it must be the case that all vertices are 
Taylor expandable to all orders in momenta [3-5]. For the purposes of this paper, 
we need only the diagrammatic rules for a particular scenario. Consider a vertex 
which is part of a complete diagram, decorated by some set of internal fields and by 
a single external A 1 (or A 2 ). The diagrammatic representation for the zeroth order 
expansion in the momentum of the external field is all that is required and is shown 
in fig. 4 [1, 16]; note the similarity to fig. 3. 




Figure 4: Diagrammatic representation of zeroth order Taylor expansion. 



The interpretation of the diagrammatics is as follows. In the first diagram on the 
right-hand side, the vertex is differentiated with respect to the momentum carried 
by the field X, whilst holding the momentum of the preceding field fixed (if the 
preceding field carries zero momentum, it is effectively transparent to the momentum 
derivative [1] and so we go in a clockwise sense to the first field which carries non- 
zero momentum to determine the momentum held constant). Of course, using our 
current diagrammatic notation, this latter field can be any of those which decorate 
the vertex, and so we sum over all possibilities. Thus, each cyclically ordered push 
forward like term has a partner, cyclically ordered pull back like term, such that the 
pair can be interpreted as 

' Vertex, (8) 



t 1 



where r and s are momenta entering the vertex. In the case that r = —s, we 
can and will drop either the push forward like term or pull back like term, since 
the combination can be expressed as 5^; we interpret the diagrammatic notation 
appropriately. 

The other diagrammatics we require is the representation of the momentum 
derivative of the effective propagators. These effective propagators are represented 
simply by a line. The momentum derivative can be with respect to the momentum 
entering / leaving either end. We indicate this by placing the momentum derivative 
symbol in the middle of the effective propagator and adding an arrow, as shown in 
fig. 5. 
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Figure 5: A differentiated effective propagator. 



Now, since the left-hand end of the effective propagator follows the momentum 
derivative in the sense indicated by the arrow, we differentiate with respect to the 
momentum entering this end. Since equal and opposite momentum enters the right- 
hand end, it is clear that we can reverse the direction of the arrow, at the expense 
of a minus sign. 

In complete diagrams, Taylor expansions of various components occur for one of 
two reasons. First, the expansion can be forced: if a diagram contains a structure 
manifestly of the order in external momentum to which we are working (in this paper, 
an 0(p 2 ) stub), then other structures carrying the external momentum can be Taylor 
expanded (so long as such a step does not generate IR divergences — see sec. 3.4.3). 
Alternatively, the expansions can be unforced, but constructed for convenience. 

In both cases, diagrams contain a discontinuity in momentum arguments. In 
the former case, the existence of an 0(p 2 ) stub makes it clear what the momentum 
routing must be. In the latter case, it is necessary to indicate where the discontinuity 
occurs. To do this, we introduce the 'bar' notation of fig. 6. 




Figure 6: Notation to indicate a discontinuity in momentum. 



Note that there is actually a choice of where we place the bar (to the left or right 
of the top vertex), depending on the momentum routing before Taylor expansion. 
However, having made a choice, the location of the bar is set, and cannot be changed. 
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2.5 Charge Conjugation Invariance 



Charge conjugation invariance can be used to simplify the diagrammatics, by allow- 
ing us to discard certain terms and to combine others. The diagrammatic prescrip- 
tion for replacing a diagram which possesses exclusively bosonic external fields with 
its charge conjugate is [1, 16] to reflect the diagram, picking up a sign for each 

1. external A 1 , 

2. performed gauge remainder, 

3. momentum derivative symbol (note that the direction of the arrow accompa- 
nying such symbols is reversed by the reflection of the diagram) . 



2.6 The Weak Coupling Expansion 
2.6.1 The Flow Equation 

Following [4, 13], the action has the weak coupling expansion 

OO _^ -i 

S = Y,{9 2 ) 1 s i = -s Q + s 1 + ---, (9) 

where So is the classical effective action and the 5i>o the ith-loop corrections. The 
seed action has a similar expansion: 

oo 

s = $> 2 *4 (10) 

Note that these definitions are consistent with S 9 = g 2 S — 25; identifying powers of 
g in the flow equation, it is clear that Si and Si will always appear together. With 
this in mind, we now define 

= S t - 2S t . (11) 
Recalling (1), the (3 functions for g and 52 are 

1 00 

Ad A - = -2^A(«)9 2(4 - 1} (12) 
1 00 

Ad A - = -2^A(l/a) ff2 2(4 - 1) , (13) 

where the Pi(ct) are determined through the renormalisation condition (5) and the 
Pi(l/a) are determined through (6). The coefficient j3\ = —f3\ is independent of a, 



13 



as we will explicitly demonstrate in sec. 3.2. For generic a, we expect the coefficient 
/32(a) to disagree with the standard value; as we will explicitly confirm in section 4.5, 
agreement is reached for /^(O). 3 

Utilising eqns. (12,13) and rewriting g\ in terms of ag\, it is apparent that Kd\a 
has the following weak coupling expansion: 



Ad\a = ^2 



2 1 



where 



7i = -2a (j3i(a)-a%(l/a)) . 



(14) 



(15) 



Substituting these definitions into (7), we obtain the weak coupling expansion of 
the flow equation, shown in fig. 7. 



{/} 



r=l L 



d 

2(n r - l)P r +lrg^ 



1 n 



r=0 




{/} 



(16) 



Figure 7: The diagrammatic form for the weak coupling flow equation. 



The symbol • = —Ad^\ a . A vertex whose argument is a letter, say n, represents 
S n . We define n r := n — r and n± := nil. The 'bar notation' of the dumbbell term 
is defined as follows: 

a o\Sn—n S r ] = CLq [S n — r , S r \ — ao[S n — r , S r ] — ao[S n — r , S r \. 
2.6.2 The Effective Propagator Relation 

The effective propagator relation [13] is central to the perturbative diagrammatic 
approach, and arises from examining the flow of all two-point, tree level vertices. 
This is done by setting n = in (16) and specialising {/} to contain two fields, as 
shown in fig. 8. We note that we can and do choose all such vertices to be single 
supertrace terms [1, 16]. 



3 We note that whilst we expect /?2(0) = /32(0), there is no reason to generically expect /32(1/q) = 
/32(a) since g and g2 are not treated symmetrically in the flow equation. 
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(17) 



Figure 8: Flow of all possible two-point, tree level vertices. 

Following [1,3-5,13,14,16], we use the freedom inherent in S by choosing the 
two-point, tree level seed action vertices equal to the corresponding Wilsonian effec- 
tive action vertices. Eqn. (17) now simplifies. Rearranging, integrating with respect 
to A and choosing the appropriate integration constants [1, 16], we arrive at the rela- 
tionship between the integrated ERG kernels — a.k.a. the effective propagators — and 
the two-point, tree level vertices shown in fig. 9. Note that we have attached the 
effective propagator, which only ever appears as an internal line, to an arbitrary 
structure. 

Af-(0) ( = M ^ _ = M ^ _ Af»^ (18) 

Figure 9: The effective propagator relation. 



The structure >> is called a gauge remainder [13]. The individual components 
of >[> will often be loosely referred to as gauge remainders; where it is necessary 
to unambiguously refer to the composite structure, we will use the terminology 'full 
gauge remainder'. 

The various components on the right-hand side of (18) can be interpreted, in the 
different sectors, according to table 1, where we assume that the gauge remainder 
carries momentum p. 



The functions f(k 2 /A 2 ) and g(k 2 / A 2 ) need never be exactly determined; rather, 
they must satisfy general constraints enforced by the requirements of proper UV reg- 
ularisation of the physical SU(N) theory and gauge invariance. We will see the effect 
induced by the latter shortly. However, we will find it useful to have concrete alge- 
braic realisations of the two-point, tree level vertices, and effective propagators (and 
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> > 


F, F 


&MN 


(M t /A 2 ,9 P ) (p„,2) 


A 1 




P^/P 2 Pv 


& 


1 





Table 1: Prescription for interpreting eqn. (18). 

hence / and g), which we collect together in appendix A. The effective propagators 
are denoted by A XY , where X and Y denote the flavours at the ends. 

It is important to note that we have defined the diagrammatics in fig. 9 such 
that there are no 1/N corrections where the effective propagator attaches to the 
vertex. We do this because, when the composite object on the left-hand side of fig. 9 
appears in actual calculations, it always occurs inside some larger diagram. It is 
straightforward to show that, in this case, the aforementioned attachment corrections 
always vanish [16]. 

2.6.3 Further Diagrammatic Identities 

We begin with a diagrammatic identity which follows from gauge invariance and the 
constraint placed on the vertices of the Wilsonian effective action by the requirement 
that the minimum of the Higgs potential is not shifted by quantum corrections: 

->(o)- =0. (19) 

This follows directly in the ^-sector, since one-point A-vertices do not exist. In 
the .F-sector, though, we are left with one-point C 1 and C 2 -vertices, but these are 
constrained to be zero. 

From the effective propagator relation and (19), two further diagrammatic iden- 
tities follow. First, consider attaching an effective propagator to the right-hand field 
in (19) and applying the effective propagator before > has acted. Diagrammatically, 
this gives 

>© = = > - »> , 

which implies the following diagrammatic identity: 

» =1. (20) 
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The effective propagator relation, together with (20), implies that 
-@ > = > - >» =0. 



In other words, the (non-zero) structure > kills a two-point, tree level vertex. 

But, by (19), this suggests that the structure > must be equal, up to some 

factor, to <\. Indeed, 

> = < , (21) 

where the dot-dash line represents the pseudo effective propagators of [1, 16]. 

The final diagrammatic identity we require follows directly from the independence 
of > on A: 

>= 0. (22) 
2.6.4 Diagrammatic Expression For (3\ 

We conclude the review by giving, in fig. 10, the diagrammatic expression of [1, 16] 
for Pi . On the right-hand side, we suppress the Lorentz indices of the external fields 
and work at 0(p 2 ), as will always be done from now on. This expression will be used 
throughout the rest of this paper, but also serves to introduce the final diagrammatic 
rules. 



4A □/«/(?) = ~\ 



Figure 10: Diagrammatic expression for Pi. 
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There are a number of things to note. First is that all diagrams are built from 
the following components: Wilsonian effective action vertices, effective propagators, 
> and >. There are no seed action vertices or covariantised kernel vertices, reflecting 
the universality of the 

Diagrams D.l and D.2 are straightforward, containing just Wilsonian effective 
action vertices and effective propagators. Diagram D.3 comprises a nested gauge 
remainder [1,16]. Diagrams D.4 and D.5 possess a common bottom vertex and 
effective propagator. This effective propagator must be in the C l -sector, else the 
diagrams vanish by charge conjugation invariance. The structures attached to the 
top end of the common effective propagator will crop up repeatedly; for convenience 
we define 

^ = ^ -4 9- . 

Diagram D.6 possesses an 0(p 2 ) stub: the renormalisation condition (5) and 
gauge invariance demand that [1, 13, 16] 

Solaip) = A p n Ma (p) = 2D Ma (p) + £>(/). (23) 

In diagrams possessing an 0(p 2 ) stub and a momentum derivative symbol, it will 
be understood that the structure attaching to the stub does not carry momentum p, 
unless explicitly indicated otherwise (see sec. 3.4.3). The structure attaching to the 
0(p 2 ) stub of dia gram D.6 is defined in terms of (pseudo) effective propagators and 
gauge remainder components: 

4 n i ^ 

: / - 2 [J ' 

Finally, we introduce some nomenclature. The complete set of diagrams inside 
the square brackets of fig. 10 will be denoted by T>i fJiU (p). The first three diagrams — 
whose sum we note is transverse [1, 16] — will be referred to as the standard set and 
the final diagram as the little set. T>\ tiv {p\ modulo the little set, will be denoted by 
v %Av)- We define 



= CO 




this expression is such that, if both external legs are A 1 s, it reduces to {$). 
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3 A-Derivative Terms 



In this section we detail the methodology used for treating A-derivative terms. Fol- 
lowing a statement of the basic idea in sec. 3.1, the methodology is developed in 
sec. 3.2 by looking at the one-loop integrals which contribute to (5\. 

In sec. 3.3, we build on the one-loop case to obtain the expected general form 
for two- loop integrals. As preparation for the evaluation of the subtraction 
techniques are introduced, which will allow us to explicitly demonstrate how non- 
universal contributions can cancel between diagrams. The methodology for this, 
which also has applications to terms which require manipulation at 0(p 2 ), is pre- 
sented in sec. 3.4. 

We conclude this section by showing that certain running couplings which can 
spoil the universality of /?2 can always be removed by a suitable choice of the seed 
action. This proves the statement made in [1] that it is possible to ensure that the 
only couplings requiring renormalisation are g and gi . 

3.1 Introduction 

The simplest A-derivative terms we will encounter are those contributing to 0\ . From 
fig. 10 we know that we can write 



We now want to make the integral over loop momentum (which we will take to be 
k) to be explicit and so write 



The next step that we wish to perform is to interchange the order of the A- 
derivative and the momentum integral. This step is trivial only if the integral is 
convergent, even after this change. Temporarily ignoring this subtlety gives: 



Since the left-hand side of this equation comprises a number times 0(p 2 ), it 
follows that the coefficient multiplying the 0(p 2 ) part of the right-hand side must 
be dimensionless. Consequently, we can schematically write 




(24) 




01 = A<9a| q (Dimensionless Quantity) . 
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For the right-hand side to survive differentiation with respect to A| a , it must ei- 
ther depend on some dimensionless running coupling — other than g and a — or there 
must be some scale, other than A, available for the construction of dimensionless 
quantities; we show how to avoid introducing such couplings in sec. 3.5. 

One scale which is available is p and so we can envisage contributions to 0\ of 
the form (in D = 4) 

Ad A |ahip 2 /A 2 . 

Indeed, the standard set gives rise to contributions of this type (see sec. 3.2.2). 
However, as we will see in sec. 3.2.3, contributions of this type cannot be formed 
from the little set — but we know from [13] that the little set does contribute to f3\. 

For the little set, then, what scales are there, other than A, available for the 
construction of dimensionless quantities? Courtesy of the SU(N\N) regularisation, 
there are no scales in the UV. Naively, we would not expect a scale to arise in the 
IR, either. However, the key point is that interchanging the order of differentiation 
with respect to Ad\\ a and loop integration in (24) has the capacity to introduce IR 
divergences. This is most easily appreciated by noting that A^ 1 ~ 1/k 2 , whereas 
A^ 1 ~ 1/A 2 (see eqn. (58)). Thus, to legally move Ad\\ a outside of the loop inte- 
gral, we must introduce some IR regulator, which then provides the scale necessary to 
form dimensionless quantities. After allowing A<9a|« to act, this unphysical scale will 
disappear. IR divergences introduced in this way will be called pseudo-divergences, 
since they are an artifact of the way in which we have chosen to perform the calcu- 
lation. 

Noting that in the case of the standard set it is effectively p which is providing 
the IR regularisation, our strategy for evaluating loop integrals is to look at the 
IR end. Scanning through the list of effective propagators (58)-(62), it is apparent 
that the leading contributions occur when all effective propagators are in the A l - 
sector; likewise for any instances of >. Gauge invariance and considerations as to 
the supertrace structure will eventually determine that all contributions to f3\ and 
(52 are ultimately limited to the lowest order momentum contributions from the A 1 
sector; it is precisely this regime — and this regime alone — that is universal. 

3.2 One-loop Integrals 
3.2.1 Vanishing Diagrams 

We begin our analysis by looking at diagrams D.4 and D.5 . Recall that these 
two diagrams have no analogue in the computation of (3\ presented in [13] (as a 

4 Although we will find a subtle interplay between the IR and UV, which is commented on at the 
end of sec. 3.2.2. See also [21]. 
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consequence of all actions being restricted to single supertrace terms) and so we had 
better find that they vanish. 

Consider the IR end of the loop integral. The bottom vertex does not carry the 
loop momentum, nor does the effective propagator leaving it. However, we know 
that this effective propagator must be in the C-sector. This immediately tells us 
that diagram D.5 is IR safe, even if we interchange the order the A-derivative and 
the momentum integral: performing this step, the loop integral just goes as 




Let us now focus on diagram D.4. To try and find IR divergences, we take 
the fields involved in the loop integral to be in the ^-sector. To deal with the 
top vertex, we recall that it is Taylor expandable in momenta [3-5]. Hence, to try 
and isolate the most IR divergent contribution from the loop integral, we take the 
minimum number of powers of momenta from the top vertex consistent with Lorentz 
invariance. Given that the field entering this vertex from beneath is in the C-sector, 
we can take O (mom ). Diagram D.4 is thus clearly IR safe in D = 4 since the loop 
integral looks at worst — without even taking into account the A-derivative — like 




in the IR. 

Hence, we can safely interchange the order of integration and differentiation with 
respect to A| a for both diagrams D.4 and D.5. Having done so, we know that (the 
0(p 2 ) parts of) both diagrams will vanish, when computed in D = 4. In preparation 
for the two-loop calculation, where the integrals are most conveniently evaluated 
using dimensional regularisation, suppose that we now work in D = 4— 2e dimensions. 
Rescaling our loop momenta k — > k/A we see that, at 0(p 2 ), the diagrams acquire 
an overall factor of A~ 2e . Working in this scheme, the A-derivative of the diagrams 
now ~ e which, of course, vanishes in the D —* 4 limit. 

Note that diagram D.l (the first element of the standard set) also vanishes in 
D = 4, after differentiation with respect to A| Q . However, we will always keep this 
term together with the other elements of the standard set. This is done because 
we often exploit the fact that the set is transverse, and want to be able to do this 
irrespective both of D and whether or not the standard set is struck by a A-derivative. 

3.2.2 IR Regularisation Provided by p 

In this section, we will encounter diagrams that survive differentiation with respect to 
A in D = 4 and for which p plays a role in the IR regularisation. The only diagrams 
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that fall into this class are diagrams D.2 and D.3. However, as just mentioned, 
diagram D.l will come along for the ride; a full understanding of the standard set is 
crucial for the two-loop calculation. 

Diagram D.l is IR safe. Diagrams D.2 and D.3, before differentiation with respect 
to A| a , have the IR structure 



where we have taken a single power of momentum from each of the three-point ver- 
tices, have chosen all effective propagators to be in the ^4 l -sector and have evaluated 
any cutoff functions at zero momentum. 

Note that choosing the effective propagators to be in the ^-sector constrains 
diagram D.2, considerably. For three-point vertices decorated exclusively by A l - 
fields, it must be the case that all fields are on the same supertrace and hence of the 
same flavour. Consequently, for the contributions to diagram D.2 with the severest 
IR behaviour, all fields must be in the ^-sector. 

Returning to eqn. (25), it is clear that the presence of p in the denominator is 
required to regularise the integral in the IR, at least when we choose to take 0(p 2 ) 
from the vertices. Performing the integral in D = 4 will then give us something of 
the form 



When this is hit by the A derivative — which we can move outside the integral — 
only the first term will survive and so we will be left with a (universal) coefficient 
multiplying two powers of p. 

That the final answer is Taylor expandable in p gives us an alternative way in 
which to evaluate A-derivative terms. Having moved the A-derivative outside the 
integral, we Taylor expand the denominator in p. Doing this, p will no longer act 
as a regulator and so we will then generate IR divergences, when we perform the 
integral. However, all divergences will be killed by the A-derivative. To parameterise 
these pseudo-divergences, we must introduce an IR regulator; it is natural to use 
dimensional regularisation. 

It may, at this stage, seem a little perverse to have traded one IR regulator, 
p — which occurs naturally — for another. However, even for diagrams which are not 
Taylor expandable in p, it will turn out that we are often interested in the Taylor 
expandable part. By Taylor expanding in p, the resulting integrals tend to be easier 
to perform. 

We now discuss how this procedure works, in more detail: 
1. Take O(mom) from each of the vertices; 




(25) 



0(p 2 )(aln(p 2 /A 2 ) + b). 
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2. Taylor expand the denominator in p; 



3. replace the upper limit of the radial integral with A, thereby cutting off modes 
above this scale; 

4. rescale k — > k/A, so that the diagram acquires an overall factor of A~ 2e . 
Having done the angular integral, we are left with an expression of the form: 

(A j f k -dkO{p 2 ). 



(4ir) D / 2 Jo k 4 

Performing the integral gives us a factor of 1/e, as expected. This is killed by a factor 
of e arising from differentiation with respect to A| a , confirming the consistency of 
the approach. 

Before moving on, we must justify the validity of the third step. We know that 
the integral we are dealing with has support only in the IR. However, J k l/k is 
not UV regularised and so we must incorporate the effects of the UV regularisation. 
Since the details of the regularisation will not affect the IR, at leading order, we 
choose the simplest form that cuts off momentum modes above the scale A. The 
non-universal corrections to this will necessarily remove any IR divergence, even 
before differentiation with respect to A| a , and so vanish in the limit that e — ► 0. 
An implicit part of this step is that we now throw away all F-sector diagrams and 
evaluate any cutoff functions at zero momentum. 

In readiness for the two-loop calculation, we will now extend our analysis of the 
standard set, and will give its general form. We know the following facts about the 
standard set: 

1. the sum of the diagrams is transverse; 

2. when struck by Ad\\ a , the coefficient of the 0{p 2 ) term is universal, up to 0(e) 
corrections (this follows from [13]); 

3. in D = 4, the diagrams have the structure 

0(p 2 )ln(p 2 /A 2 ) + 0(p 2 ) + ---, 
where the ellipsis denotes terms which are higher order in p. 

From these three points and dimensions it follows that, in D = 4 — 2e, the standard 
set takes the algebraic form 



NA 



-2c 



(4tt) d / 2 
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n a p(p) + 0(p\p 
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where ao is a universal coefficient and the ellipsis denotes terms which are higher 
order in e. 

We have no further interest in the 0(p ,p 4 ~ 2e ) terms and so turn to the ellipsis. 
There are two ways in which we expect these terms to arise. On the one hand, 
compared to eqn. (25), we can take additional powers of k in the numerator of the 
integrand, giving us non-universal contributions which are Taylor expandable in p. 
On the other hand, we generically expect the coefficient ao to have arisen from some 
function of D in which we have taken D = 4. Expanding this function in e will give 
rise to sub-leading contributions. These contributions will be called computable. At 
the two-loop level, we will see how computable contributions can combine to give 
universal quantities. 

Just as we talk of computable parts of some diagram, so too will we talk of 
the complimentary non-computable parts. We emphasise that by non-computable 
we really mean that the corresponding coefficients cannot be computed without 
specifying non-universal details of our set-up (i.e. cut-off functions and the precise 
form of the covariantisation) ; it is not that we cannot, in principle, calculate them. 

Hence, the standard set takes the following form 

O a p(p) + 0(p\p 4 - 2 *), (26) 

where bo = — ao, the aj>o are a mixture of computable and non-computable contri- 
butions and the bi are entirely computable. 

Notice that one-loop computations are insensitive to the p~ 2e /A~ 2e terms: taking 
into account the additional factor of A _2<E sitting outside, such terms are independent 
of A and so will be killed by the A-derivative. At two-loops, where the standard set 
can occur as a sub-diagram, we expect such contributions to survive. 

We now discuss how to calculate the coefficients bi. It is convenient to begin by 
contracting the standard set with 8 a p. The 6, arise from integrals of the form 

f 0(p 2 ,p.k,k 2 ) 

Jk k 2 (k-p) 2 • 

Note that the 0(k 2 ) term does not contribute to the bi\ the denominator becomes 
just l/(k—p) 2 and so, by shifting momentum, we can remove p from the denominator 
entirely. It is now not possible to generate a power of p~ 2e . 

The next step is to combine denominators, using Feynman parameterisation: 

I dah{a) J k (k 2 + K 2 r 



-2c 
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where K 2 = a(l — a)p 2 and h(a) is some function of a. At this stage, it is now tempt- 
ing to proceed as before and restrict the range of the radial integral. However, this 
leaves us with an unpleasant calculation as we cannot use standard dimensional reg- 
ularisation formulae. Besides, there is a much simpler way to proceed: differentiate 
twice, with respect to p 2 . 

The effect on the integral is to ensure that it is UV regularised by the denomi- 
nator of the integrand — without the need for any cutoff regularisation. We call this 
automatic UV regularisation (which will play an important role at two loops). Since 
we are interested in the part of the integral which has support in the IR, there is 
no need for us to restrict the range of integration, as doing so would only serve to 
make the calculation harder. Retaining just A-sector diagrams, we evaluate all cutoff 
functions at zero momentum, leaving us with an integral we can do using standard 
dimensional regularisation techniques. 

We perform the integral and compare it to the second derivative with respect to 
p 2 of eqn. (26), contracted with 5 a p. This is one place where the value of keeping 
the standard set together manifests itself: because we know the standard set to 
be transverse, we know the effect of contracting with 6 a g. Equating powers of e 
allows us to determine the 6,. The first two coefficients, computed in dimensional 
regularisation are: 



where 7em is the Euler-Mascheroni constant. Note that &o is n °t just computable 
but also universal, being independent of the way in which we choose to evaluate 



At this point it is worth pausing to consider why the coefficients ao and 6, have no 
dependence on N. We note in eqn. (26) that we have extracted an overall factor of N, 
but we might suspect that the a, and bi incorporate attachment corrections. 5 Let us 
look first at diagram D.3. Since this is formed by the action of gauge remainders on 
three-point (tree level) vertices decorated by an external field, we know from [1, 16] 
that we can discard all 1/N corrections. 

Now consider diagram D.2. The contributions to ao and &, come when all fields 
are in the ^-sector. If either of the effective propagators attaches via a 1/N correc- 
tion, then the external fields are always guaranteed to be on the same supertrace, 
irrespective of location: the diagram vanishes by charge conjugation, since a three- 
point A 1 vertex changes sign under the interchange of any pair of fields. However, 

5 The overall factor of N is what we expect for the diagrams if there are no attachment corrections. 
In this case, each diagram comprises one loop decorated by the two external fields and one empty 
loop. The empty loop yields str<r+ = N if the internal fields are bosonic and strcr_ = —N if the 
internal fields are fermionic. 



b 
bi 
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124/9 + 20 7EM /3 



(27) 
(28) 



it [13]. 
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we expect the etj>i, to be non-trivial functions of 1/N, since these coefficients receive 
contributions in which the vertices of diagram D.2 are each charge conjugation even 
(i.e. AAC vertices) and from diagram D.f. 

We conclude our analysis of the standard set by noting a beautiful interplay 
between the IR and the UV, illustrated by diagram D.3. The strategy we have used 
is to pull the A-derivative outside of the momentum integral and then focus on the 
IR end. Focusing on the IR end allows us to throw away the regularising diagram. 
However, the regularising diagram was required to define the ^-sector diagram when 
we interchanged the order of differentiation with respect to A and loop integration. 

Now, suppose that we had left the A-derivative inside the integral. Then the 
^-sector diagram actually dies, since ^-sector gauge remainders are independent 
of A. We are left with the -B-sector diagram, which provides the same leading order 
contribution as the ^-sector diagram, but arising from the UV! 

Interplay such as this will only arise when the components of some diagram which 
gives a contribution in the IR are not regularised by cutoff functions, alone (c/. [21]). 

3.2.3 Loop Integrals Independent of p 

We conclude our survey of one-loop integrals by looking at the final diagram which 
contributes to f3i, D.6. The difference between this diagram and the ones just anal- 
ysed is that p is not involved in the IR regularisation. The integral just goes like 

Ik 

in the IR and we can use the techniques of the previous section to evaluate such 
terms. We thus expect the complete set of diagrams contributing to 0\ to take the 
following form, before differentiation with respect to A| a : 

a^p) + 0(p\p^), (29) 

where the Bi are computable and the Ai generally contain both computable and 
non-computable parts. The universal coefficient Aq yields the sole contribution to 
j3i, in the e — > limit. 

3.3 Two-loop Integrals 

In this section, we develop the machinery of the previous section to deal with two- 
loop A-derivative terms. The integrals we have to deal with fall into two classes: 
factorisable and non-factorisable. In the former case — dealt with in sec. 3.3.1 — the 
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loop-integrals are independent, whereas, in the latter case — dealt with in sec. 3.3.2 — 
they are not. 

Following on from the one-loop case we expect and, indeed, find [16] that we can 
write 6 

4&n^(p) = -£[xv(p)]"- 

One of the main sources of complication in the two- loop case is that, even after 
differentiation with respect to Ad\\ a , individual elements of V^^vip) can still possess 
IR divergences. It is only the sum of diagrams contributing to ^^{p) that we expect 
to give a finite (universal) contribution after differentiation with respect to Ac^j^. 

Since we are interested in two-loop integrals which contribute to we will, up 
to factors of p~ 2e , work at 0(p 2 ). 



3.3.1 The Factorisable Case 

To understand the algebraic form of two-loop diagrams, we need first to understand 
their structure. Since we are dealing with factorisable terms, we expect them to 
comprise two one- loop sub-diagrams, each of which carries external momentum p. 
These sub-diagrams must be connected to each other, and so we predict that they 
will be joined together by an effective propagator (for explicit examples, see sec. 4.1). 
This effective propagator just contributes powers of the external momentum and so 
we take the general form of a factorisable two-loop integral to be: 

where we obtain a power of p~ 2e /A~ 2e for each loop in which p provides IR regular- 
isation. 



3.3.2 The Non-Factorisable Case 

We expect non-factorisable diagrams to possess the same ingredients as factorisable 
ones, but joined together in a different way. Taking the loop momenta to be I and 
k, we know that one of the effective propagators must ~ 1/(1 — k) 2 (assuming it to 
be in the ^-sector). Conservation of four- momentum at a vertex then implies that 
there must be at least one effective propagator carrying I and at least one carrying 
k. Knowing that, at 0(p 2 ) the integrand must be of mass dimension —8, we expect 
the most divergent type of diagram we can construct (assuming no IR regularisation 

6 In the limit that a — > 0. 



27 



is provided by p) to take the following form in the IR: 




0(mom 2 ) 



i )fc k 2 {l-k)H^ 



In fact, we will see in sec. 4.4 that, for the set of diagrams contributing to 
gauge invariance prevents the appearance of 1/Z 6 and so we would find that, taking 
the above form, we would be forced to have I 2 in the numerator. Hence, the most 
divergent type of integral we find has the following structure in the IR: 



To evaluate the contribution coming from the IR, we observe that the /-integral is 
automatically UV regularised. Thus, using dimensional regularisation, we perform 
the Z-integral first, with unrestricted range of integration, and perform the /c-integral 
second, with the range of radial integration restricted to A. We obtain one power 
of 1 /e from the Feynman parameter integral and a further power from the radial k- 
integral. Doing the integrals the other way around would be awkward, as we cannot 
then use an un-restricted range of integration for the inner integral. 

Had there been ps present in the denominators, providing regularisation, we 
would expect accompanying factors of p~ 2e /A~ 2e . Consequently, eqn. (30) is the 
form for a generic two-loop integral. 

3.3.3 Considerations for fa 

For the actual computation of fo, we can constrain some of the coefficients in 
eqn. (30). Since the terms corresponding to the coefficients ej are independent of 
A, we can drop them as they will vanish after differentiation. Given that we are 
comparing our final answer with the Taylor expandable, finite expression /^□^(p), 
it must be that the coefficients cq and di vanish. The coefficients cq and do are 
entirely computable. The coefficient d\, on the other hand, comprises both a purely 
computable part and a non-computable part multiplied by a computable coefficient. 
The computable part and the computable coefficient must both be zero. 

Ultimately, we will be left with the coefficient c\ being the only contribution to 
the final answer. One of the primary tasks ahead is to show that the non-universal 
contributions to c\ cancel between diagrams. This problem has two sides. First, we 
must show that non-computable contributions from vertices etc. cancel out. Then 
we must show that the computable contributions to c\ combine to give the standard, 
universal answer. 




(31) 
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3.4 Subtraction Techniques 



3.4.1 Basics 

Rather than attempting to process two-loop diagrams directly, we perform an in- 
termediate step whereby we add and subtract a set of terms designed to remove 
all non-computable contributions from the calculation. We illustrate this technique 
with a simple example. Consider the two-loop integral arising from the computation 
of the scalar two-loop /^-function, within the ERG of [14]. 



where we call the second term a subtraction and the third term its corresponding 
addition. The addition trivially combines with the final term, though this is specific 
to this example and of no particular significance. In deference to the forthcoming 
gauge theory calculation, we will leave the third and fourth terms uncombined, so 
that we can illustrate the strategy generally employed. 

Let us start by focusing on the k integral, in the first term. Following [14,22], 
we know that we can set c(l — k) = c(k), as contributions higher order in I are 
killed in the e — ► limit. Next, use the by now familiar prescription for the cutoff 
functions: if the integral is regularised without the cutoff functions, then we simply 
evaluate them at zero momentum, leaving the domain of integration unrestricted. If 
the integral requires the cutoff functions for regularisation, then restrict the domain 
of integration and evaluate the cutoff function at zero momentum. The sub-leading 
corrections to this will manifest themselves as additional powers of momenta, in the 
numerator. 

If we were to start taking such sub-leading (non-computable) contributions under 
the fe-integral, then this will allow us to Taylor expand the k- integral in I. For the 
Z-integral to still diverge in the IR — and thus survive differentiation with respect to 
A| Q (in the limit that e — » 0) — we must take 1°. However, such non-computable 
contributions will be cancelled by exactly the same contributions coming from the 
second term, above. Hence we have: 



where C tells us that, when considered as a pair, the first two fc-integrals yield a 
computable contribution. 




(32) 



where A/ = c(l)/l 



2 . We can trivially rewrite this as 
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For the first term not to die in the e — > limit, we must set c(l) — > 1 and so 
can extend the 'C to cover the whole term. The second term, derived from the 
subtraction, can now be combined with the addition. Generally, we will always 
perform this step: it isolates non-computable contributions of the first term in (32) 
that survive in the D — > 4 limit. We thus have: 



l,k 



A^_ 



A! 



1 



NC 



A?A2 



+ 0(e) 



(33) 



The final step is to re-express each of the components of the last term as C + NC 
and multiply out, collecting terms by using the freedom to interchange I and k. 



Noting that the 
purely computable contribution: 



(a?; 



NC 



NC 



term vanishes in the D — > 4 limit we are left with a 



+ 0{e) 



We have thus demonstrated that the original integral does, indeed, give something 
which is computable. 

To calculate this integral, we can use a mixture of the techniques already dis- 
cussed. In the factorisable case, we simply restrict the ranges of the integrals. In the 
non-factorisable case, we note that the /-integral is automatically UV regularised. 
Hence, we perform this integral first with an unrestricted domain of integration but 
then restrict the domain for the remaining A;-integral. It is straightforward to confirm 
that we reproduce the expected, universal answer 
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3.4.2 Generalisation to the Gauge Case 

Constructing subtractions in the gauge case is exactly analogous to the simpler case 
just analysed, but the structure of the cancellations of non-computable contributions 
is much richer. In the same way, the subtractions are constructed such as to isolate 
non-computable contributions, by noting that denominators can be Taylor expanded 
in momenta if sufficient powers of momenta are present in the numerator. 

As the whole formalism is based around Taylor expansion, it is not surprising 
that we will need to use the techniques of sec. 2.4 to Taylor expand vertices in 
momenta. We know that the lowest order terms constitute momentum derivatives 
of lower point vertices, and that the sign of this derivative depends on whether we 
have had to push forward or pull back. We define the subtraction to be the term 
which removes non-computable components from the parent diagram, and not by its 
sign. Hence, a subtraction involving a pull back will come with a positive sign. 
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The real subtleties in the gauge case arise because the procedure of constructing 
subtractions and additions and the subsequent cancellation of all non-computable 
contributions leads, generically, to a complete loss of momentum routing invariance. 
As we will see, this invariance can be partially restored, though it is not necessarily 
efficient to do so. 

As a first illustration, suppose that we have taken a non-factorisable, two-loop 
diagram and constructed a factorisable subtraction. Putting the addition to one side 
for the moment, we will suppose that the sub-diagram of this latter term, to which 
we apply C, is just the diagram shown in fig. 11. 



D.7 




Figure 11: A sub-diagram, to which we will apply C, of some factorisable subtraction. 



The algebraic form of diagram D.7 is 



□ M/3 (p)(A- 2 T / d k a 



k 2 



k 2 ' 



(34) 



where = (c/. eqn. (23)). To compute the part of this diagram left over, after 
combining with our non-factorisable term we do the following: evaluate the cutoff 
function at zero momentum, restrict the range of the integral, and proceed as usual. 

Next, suppose that we were to move the momentum derivative from the B^/k 2 
term to the kp/k 2 term, throwing away the total derivative, in the process. This 
yields 



-□tf(p)(A-*)-||f# 



kp_ 
k 2 



(35) 



Evaluating this integral in the usual way gives a different contribution, at sub- 
leading order, than the integral of eqn. (34). What is going on? The point is, 
that whilst going from eqn. (34) to eqn. (35) is usually a perfectly valid step, it 
breaks down when these terms are under the influence of C. Specifically, because the 
effect of C has been to replace cutoff functions with a restricted range of integration, 
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we are no longer justified in throwing away what would previously have been total 
momentum derivative terms. Equivalently, we have lost the freedom to shift k, under 
the integral. 

If we wish, momentum routing invariance can be retained, in this case, by modi- 
fying what we mean by C. We can use the following prescription: reinstate a term to 
both the parent and subtraction (after they are under the influence of C) such that, 
at sub-leading order, we can move momentum derivatives around with impunity. To 
understand what this term must be, let us return to eqn. (34). Rather than discard- 
ing the cutoff function straight away, we will first allow the momentum derivative to 
act. 

Doing so, averaging over angle and substituting x = k 2 yields: 



where ftp is the area of the unit sphere is D dimensions divided by (2tv) d . Integrat- 
ing the first term by parts, and discarding the resulting surface term gives: 



Now if we remove the cutoff function and restrict the range of integration, it is 
apparent that the B' term has provided a sub-leading contribution; indeed, this is 
precisely the sub-leading contribution we are after! 

We have seen how, when cutoff functions are necessary for UV regularisation, 
their derivatives can supply sub-leading contributions in the IR. This can be rephrased 
by saying that, under the influence of C, total momentum derivative contributions 
can no longer be discarded, at sub-leading order, unless we reinstate terms to parent 
and subtraction. 

When dealing with automatically regularised integrals, however, total momen- 
tum derivatives can be thrown away, even under the influence of C. This follows 
because the range of integration need not be restricted, even after we have evaluated 
any cutoff functions at zero momentum. Equivalently, in this case, we can move 
momentum derivatives around without the need to reinstate the derivative of cutoff 
functions. 

Had we combined the parent, addition and subtraction as in the scalar exam- 
ple, similar considerations to those above would apply, as follows directly from the 
complementarity of C and NC. 

Unfortunately, reinstating cutoff functions in this way can create more problems 
than it solves. To understand why, let us consider a second example, which we will 
encounter in the evaluation of Our starting point is the diagram on the left-hand 
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side of fig. 12. On the right-hand side is an exactly equivalent representation: we 
have simply taken two copies of the diagram and reflected one of them using charge 
conjugation invariance [1, 16]. In readiness for the construction of the subtractions, 
we have chosen differing momentum routings for the two diagrams. (Note, though, 
that at this stage the two diagrams are exactly equivalent, since complete momentum 
routing invariance has not been broken.) 



D.8 D.9 




Figure 12: A diagram for which subtractions will be constructed. 

As a consequence of our choices of momentum routing, diagrams D.8 and D.9 
will have different subtractions. We do this for convenience, as the subset of terms 
thus generated can be manifestly combined into something useful. The subtractions 
are shown in fig. 13, where we have used charge conjugation to collect terms; each 
diagram possesses two labels: the top one for the subtraction and the bottom one 
for the corresponding addition. 

Note that, in diagrams D.14, D.20 and D.22 (and the corresponding additions), 
the field attached to the circle representing the momentum derivative is implicitly in 
the J 4 1 -sector. This follows because the only derivatives of vertices we consider are 
those which arise from a field in the A % -sector carrying zero momentum (and here, 
the A 1 can only be an A 1 ). 

To begin the analysis of diagram D.8 and its subtractions (D.10, D.12 and D.14), 
we focus on the three-legged sub-diagram carrying loop momentum k (the momen- 
tum routing of subtractions follows from the parent). For the diagram as a whole 
to have any chance of surviving in the e — > limit, the (internal) legs leaving the 
sub-diagram must be in the ^-sector. Thus, by Lorentz invariance, the sub-diagram 
carrying loop momentum k must go as odd powers of momentum (up to additional 
non- Taylor expandable functions of I). Noting that, in D = 4, the sub-diagram car- 
rying loop momentum k goes as, at worst, (In I) x 0(mom, . . .) in the IR, it is clear 
that we must take only the O(mom) part of the sub-diagram. 
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D.14 
D.15 




D.20 
D.21 




D.22 
D.23 



T2 




Figure 13: Subtractions for diagrams D.8 and D.9. 



The effect of diagrams D.10 and D.12 is now immediately clear: they completely 
remove from diagram D.8 all contributions in which the sub-diagram carrying loop 
momentum k goes as I. 

Let us now suppose that we take 0{p) from the sub-diagram carrying loop mo- 
mentum k. We start by noting that the only place for this power of p to come from 
is the four-point vertex. Now, if all the fields leaving the four-point vertex are in the 
^-sector, then Lorentz invariance forces us to take an additional power of momen- 
tum from the four-point vertex. Recalling that we should not take any further powers 
of / or p, we see that we must take (at least) one power of k from the four-point 
vertex (and a further power of k from the three-point vertex) . The /c-integral is now 
Taylor expandable in /. In the case that the fc-dependent fields leaving the four-point 
vertex are not in the ^-sector, the /c-integral is trivially Taylor expandable in I. 

The remaining subtraction, diagram D.14, removes all these contributions; hence 
diagram D.8 turns out to be completely cancelled by its subtractions, (up to contri- 
butions which die in the D — > 4 limit). 

The same cannot, however, be said for diagram D.9. It is clear that whilst 
diagrams D.16 and D.18 remove from D.9 all contributions from the /c-integral that 
are Taylor expandable in p, these diagrams possess components which are not Taylor 
expandable mp and so will survive. Diagrams D.20 and D.22 remove all contributions 
that are Taylor expandable in I (as always, this statement is correct only up to 
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contributions that vanish anyway in the e — > limit). Non-computable contributions 
from the £>integral are precisely those which are Taylor expandable in I and p and 
so cancel between the parent diagram and its subtractions. 

That diagram D.8 is completely cancelled by its subtractions is in some senses 
a coincidence; it is certainly not generically true that parent diagrams are thus 
cancelled. Indeed, it is most efficient not to make use of this fact. Thus, to proceed, 
we employ the strategy outlined earlier: first combine parent with subtraction and 
then combine the remaining components of the subtraction with the corresponding 
addition. This yields the diagrams of fig. 14, up to 0(e) corrections. 

There are a number of important comments to make about the diagrams of 
fig. 14. First, great care must be taken interpreting the precise meaning of C and 
NC. Consider evaluating diagram D.24, for which all fields must be in the A 1 sector. 
It is most convenient to make the /-integral the inner one, as this is automatically UV 
regularised. Now, since we discard all cutoff functions (and regularising diagrams), 
the ^-integral is no longer invariant under shifts of k. This is crucially important. 

Notice that the four-point vertex attaches to the three-point vertex on its left 
via two effective propagators carrying k and I — k. We suppose that the dummy 
indices associated with these effective propagators are p and a, respectively, and 
that the other index of the three-point vertex is a. Recall that the diagrammatic 
representation used for the three-point vertex actually stands for all independent, 
cyclically ordered, three-point vertices. Since three-point A 1 vertices only exist as 
single supertrace terms, there are two independent forms for our three-point vertex: 

Saap k — I, —k) and S^p^(l, —k, k — I). 

Naively, it should be the case that the value of the diagram is the same in both 
cases; after all, going from one cyclically ordered three-point to the other is equiv- 
alent to relabelling the (internal) effective propagators. However, this relabelling is 
achieved via the shift in momentum k — > —k + l and so we know that the value of the 
diagram, under the influence of C, will actually be different for each of the cyclically 
ordered three-point vertices. This explains why diagrams D.26 and D.27 have been 
kept separate, despite the fact that these diagrams could be directly combined if 
they were not under the influence of NC (indeed, precisely such a combination was 
formed when constructing the subtraction D.12). 

We could remove the subtleties associated with a lack of momentum rerouting 
invariance, order by order in e, by re-instating the necessary UV regularisation in 
the k integral for every diagram of fig. 14. This amounts to a different prescription 
for what we mean by C and will generally yield different answers for computable 
diagrams. Of course, it is not possible to unambiguously extract the computable 
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component of a non-universal diagram; only a universal quantity formed from a sum 
of computable components is independent of the computational scheme. However, 
making this change to the prescription for C is not generally a good idea, as it 
interferes with a crucial diagrammatic step. Suppose for a moment that the top- 
most sub-diagram of D.28 were not under the influence of NC. Then this diagram 
could be immediately processed using the effective propagator relation: 




Can this manipulation be performed under the influence of NC? To answer this, it 
is simplest to think of NC acting on a diagram as the full diagram, minus C. Clearly, 
the full diagram can be processed. Now, the effective propagator relation holds if 
we take the lowest momentum contribution of both the (^-sector) two-point, tree 
level vertex and the corresponding effective propagator i.e. the relation works under 
C if all cutoff functions have been thrown away. Though by no means the only way 
to proceed, it turns out that the most efficient is to completely sacrifice momentum 
routing invariance, thereby allowing the effective propagator relation to be straight- 
forwardly applied. The processing of gauge remainders is unaffected by C and hence 
NC. The only diagrammatic step that cannot be made under the influence of our 
current prescription for C, NC is that of throwing away total momentum deriva- 
tive terms sitting under outer integrals. Note, though, that such terms can always 
be reduced to computable terms (at two loops). Consider a diagram containing a 
sub-diagram both struck by a momentum derivative and under the influence of NC. 
Since a total momentum derivative of a diagram not under the influence of either 
C or NC vanishes, it follows that, in this case, the NC can be replaced by C, at the 
expense of a minus sign. Now, since such total momentum derivative terms only 
contribute at sub-leading order in e then, in order for the diagram as a whole to 
survive differentiation with respect to A\ a in the e — > limit, we must extend the 
effect of C to cover the entire diagram. 

The second point to make about the diagrams of fig. 14 is that the sub-diagrams 
of diagrams D.25-D.27 and D.30-D.32 carrying loop momentum k can very nearly 
be combined into momentum derivatives of the second element of the standard set. 7 
Had we not originally split the parent diagram up, as in fig. 12, this fact — which turns 
out to be extremely useful in the evaluation of — would not have been manifestly 
obvious. 

Let us examine diagrams D.25-D.27 more closely (an identical analysis can be 
performed for D.30-D.32). From sec. 2.4, the component of diagram D.25 in which 

7 The external fields of these sub-diagrams must be in the J 4 1 -sector for the diagrams as a whole 
to survive in the D — > 4 limit. 
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the cyclically ordered momentum arguments of the top-right vertex are (k, I — k, —I) 
combines with diagram D.26 such that, if the overall sign of the diagram is taken to 
be positive, then the top-right vertex is struck by 



-d~ 



(37) 



Similarly, the component of diagram D.25 in which the cyclically ordered momentum 
arguments of the top-right vertex are (I — k, k, —I) combines with diagram D.26 such 
that the top-right vertex is struck by 



l-k 







l-k 







(38) 



Eqns. (37,38) differ by the term d l 



If diagrams D.25-D.27 were not under the 
influence of NC, this apparent difference would vanish, as a consequence of momen- 
tum rerouting invariance. In the current case, it survives, of course, though it can 
be converted into an entirely computable term, at the expense of the usual minus 
sign. 

Diagrams D.25-D.27 and D.30-D.32 naturally combine with the term shown in 
fig. 15, which is generated by the manipulation of a different parent diagram from 
that of fig. 12 (this parent turns out to be diagram T.13, which we will encounter in 
sec. 4.1, when we examine the diagrams contributing to fo)- 



D.36 



2 ^J?^ 




l-p 



NC 



Figure 15: The partner term to diagrams D.25-D.27 and D.30-D.32. 



If it were the case that diagrams D.25-D.27, D.30-D.32 and D.36 were not under 
the influence of NC, then they could be combined into a single term, in which the 
fc-dependent sub-diagrams sum up to give d l u ^ of the second element of the standard 

set. (Note that if the d l u ^ strikes the effective propagator carrying I — k, then it 

can be directly exchanged for — , since the effective propagator depends only on 
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(I — k) 2 . Such a term vanishes as a consequence of momentum rerouting invariance.) 
We represent the derivative of the second element of the standard set and, more 
generally, T>^ , with respect to its external momentum by 



In the case that PfL, is under the influence of either C or NC, the interpretation 
of this notation changes slightly for the second and third elements of the standard set 
(recall that all other elements of the standard set are entirely non-computable). For 
the second element of the standard set, the interpretation is taken to be consistent 
with diagrams D.25-D.27, D.30-D.32 and D.36 i.e. eqn. (39) now includes derivatives 
of the vertices with respect to the internal momentum of the diagram, according 
to (37,38). For the momentum derivative of the third element of the standard set 
we take: 



Were this pair not under the influence of either C or NC, then they could be combined. 
3.4.3 Application to Terms Manipulable at 0(p 2 ) 

The application of the subtraction techniques we have described is not limited to 
A-derivative terms, but is useful for any class of terms in which we wish to perform 
Taylor expansions. We have already encountered such manipulations when we dealt 
with the 0(p 2 ) terms in the f3\ diagrammatics of [1, 16]. For example, the elements 
of the little set were derived from diagrams which were Taylor expanded in p. Let us 
return to this by considering one of the parents of the little set, which is reproduced 
in fig. 16. 

Due to the fact that the structure carrying momentum I — p is an undecorated 
kernel, rather than an effective propagator, it is clear that p is not required to 
regularise this diagram in the IR; hence we can expand not only the vertex but also 
the kernel to zeroth order in p. For this diagram, there is no need to construct a 
subtraction. Indeed, at one-loop, there is never any need to construct subtractions 
for diagrams manipulable at 0(p 2 ). 

At two- loops, however, the situation is different. Consider the first diagram 
shown in fig. 17; in anticipation of what follows, we have constructed a subtraction. 




(39) 
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Figure 16: A Taylor expandable, one-loop diagram. 
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D.40 
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Figure 17: A two-loop diagram with an 0(p 2 ) stub, which cannot be Taylor expanded 
in p, and its subtraction. 

We begin by focusing on diagram D.38. Since we can always Taylor expand 
vertices in momenta, let us suppose that we take a power of I from the top-most 
vertex (we cannot take any powers of p, at 0{p 2 )) and let us choose to take a power 
of k from the other vertex. The leading IR behaviour of the /-integral is now 



this is not Taylor expandable in p. Note that had we taken a power of / from the 
right-hand vertex, rather than a power of k, then the extra power of I in the integrand 
would render the diagram Taylor expandable in p (to the required order). 

Now let us consider the subtraction and addition. The addition (diagram D.40) 
is manipulated in the usual way; this is basically what we would like to have done 
with diagram D.38, in the first place. The effect of the subtraction on the parent is 
to cancel all those components which are Taylor expandable in p. This immediately 
tells us the following about any surviving contributions to diagram D.38: 
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1. all fields carrying momentum I must be in the A -sector; 

2. we must take 0(1°) from the /c-integral (note that the ^-integral is Taylor 
expandable in I); 

3. we must discard any remaining contributions to the I integral which do not 
The contributions to diagram D.38 not removed by diagram D.39 are shown in fig. 18. 



D.41 




Figure 18: The contribution to diagram D.38 not removed by its subtraction. 

As required, we have taken the O(l ) from the /c-integral. The tag p~ 2e demands 
that we take the p~ 2e component of the diagram. Note, of course, that this tag 
implicitly assumes that we are using dimensional regularisation. However, were we 
to use some other means of regularising IR divergences, diagrams such as D.41 would 
still exist, but the tag would be appropriately generalised. 

3.5 Ensuring Universality 

The central tenet of our analysis of the A-derivative terms has been that the deriva- 
tive with respect to A| a of a dimensionless integral must vanish, unless there is a 
scale other than A, from which we can construct dimensionless quantities. Implicit in 
this is that there are no dimensionless running 8 couplings, hidden in the integrand. 

The most obvious candidates for dimensionless running couplings can immedi- 
ately be discounted: g counts the loop order, and so never appears in loop integrals 
and the presence of a is irrelevant, since it is held constant, when differentiating 
with respect to A. The first question we must address is whether there are actually 
any other candidates for dimensionless running couplings. 

8 Where, strictly, we mean running with respect to A| a . 
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To see how they could arise, in principle, consider the flow of any vertex, with 
mass dimension > 0. Now Taylor expand in momenta and focus on the term which 
is the same order in momenta as the mass dimension of the vertex. The coefficient 
of this term must be dimensionless; if this coefficient flows, then we have found what 
we are looking for. 

As a first example, let us consider the flow of an m-loop vertex, decorated by an 
arbitrary number, q, of C l s. We take the C l s to carry momenta rj. Recalling that 
C*s are of mass dimension zero [5, 13], all such vertices are of mass dimension four. 
Hence, we are interested in the O (mom 4 ) component of each of the vertices. 

The crucial point for what follows is that, no matter what the value of m, the 
flow is guaranteed to produce a certain type of term: specifically, we will always 
have a dumbbell structure consisting of a two-point, tree level vertex, joined by an 
undecorated kernel to a seed action vertex. This is illustrated in fig. 19. 



D.42 D.43 

=- X 

r> (0) 
! t n 

Figure 19: The flow of a vertex decorated by an arbitrary number of C l s. 

The first ellipsis denotes diagrams of the same structure as D.42 and D.43 but 
for which a different C l decorates the two-point, tree level vertex. Each of these 
diagrams possess a seed action vertex. These seed action vertices are the highest 
loop vertices which appear; moreover, all other vertices generated at this loop order 
possess fewer legs. The second ellipsis denotes the remaining terms generated by the 
flow. 

Focusing on the O (mom 4 ) components of all diagrams generated by the flow, we 
now tune the m-loop, (/-point, seed action vertices to exactly cancel the remaining 
terms. This choice of seed action is one we are entirely at liberty to make; it ensures 
that there are no hidden running couplings in this sector of the calculation. It perhaps 
seems a little artificial that we only ensure universality after some (implicit) choice 
for the seed action. We must remember, though, that /3-function coefficients are not 
strictly universal, and that scheme dependence even at one-loop is not necessarily 
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a sign of a sick formalism [13]. Our choice of seed action is merely done to allow 
comparison of the values we compute for f3\ and P2 with those computed in, say, 
MS. 

In anticipation of what follows, we emphasise that the crucial ingredient in what 
we have just done is that the flow of an m-loop, q-point vertex generates an m-loop, 
(/-point seed action vertex. Moreover, there are no vertices generated with higher 
loop order, and for the rest of the same loop order, the number of legs is < q. 
Consequently, for each Wilsonian effective action vertex whose flow we compute, it 
is a different seed action vertex we tune. This ensures that we are never in the 
situation where we have to try and tune the same seed action vertex in two different 
directions. 

Let us now move on to consider an m-loop vertex decorated by q C l s and also 
by a single A . The vertex is now of mass dimension three and so it is the O (mom 3 ) 
part we are interested in. This time, we are guaranteed to generate m-loop, q + 1- 
point seed action vertices, joined to a two-point, tree-level vertex by an undecorated 
kernel. Now, however, we see a difference between this case and the previous one: the 
tree- level, two-point vertex can be decorated by either an A 1 or a C". We illustrate 
this in fig. 20, where we take the A 1 to carry momentum p and the C l s to carry 
momenta rj. 



D.44 D.45 




Figure 20: The flow of a vertex decorated by a single A and an arbitrary number 
of 0*8. 

The first ellipsis denotes diagrams of the same structure as D.44 but for which a 
different C l decorates the two-point, tree level vertex. The final ellipsis denotes the 
remaining terms generated by the flow equation. 

Having decorated with an A 1 , we must now take account of gauge invariance. 
The most obvious effect is that the two-point, tree level vertex of diagram D.45, 
unlike that of diagram D.44, is forced to be 0(p 2 ). Given that we are working at 
O (mom 3 ) , this means that we must take a single power of momentum from the 
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corresponding m-loop, seed action vertex. Now, if this single power of momentum is 
p then, by Lorentz invariance, its index must be contracted with the two-point, tree 
level vertex — killing it. Hence, this single power of momentum must be one of the 
r.j. In turn, this means that we can Taylor expand the m-loop, seed action vertex of 
diagram D.45 to zeroth order in p; the effect of this is to reduce it to the derivative of 
a (/-point vertex. This means that gauge invariance has caused us to lose one of our 
m-loop, q + 1-point seed action vertices. We can still tune the flow of our Wilsonian 
effective action vertex to zero, though, by virtue of the presence of diagram D.44 
and the diagrams represented by the first ellipsis. 

Let us examine the tuning of the seed action vertices in a little more detail. To do 
this, consider contracting the diagrams of fig. 20 with the momentum carried by the 
A 1 . On the left-hand side of the equation, we now have the flow of (a set of) m-loop 
vertices, decorated by q C l s. Since we were working at O (mom 3 ) but have contracted 
with a power of momentum, we should now be looking at O (mom 4 ) . However, we 
know from our work on pure-C l vertices that the flow of such terms has already been 
tuned to zero. Thus, returning to the diagrams of fig. 20, gauge invariance ensures 
that we need only tune the seed action vertices so as to remove those O (mom 3 ) 
contributions transverse in p\ Note that this reproduces the conclusions of [13], in 
which the special case of the flow of a tree level ACC vertex was considered. 

Next, we extend our analysis to a vertex decorated by two A l s and q C l s. We 
now work at O (mom 2 ). If we take q > 0, then our analysis just mirrors what we 
have done: to avoid dimensionless, running couplings, gauge invariance ensures that 
we need only tune the seed action vertices transverse in the momenta of the A l s. 
What if q=0? Now the renormalisation condition guarantees that the flow of the 
vertex vanishes; this is, of course, exactly what we have been utilising to compute 
/3-function coefficients. 

In the case of a vertex decorated by three ^4*s and q C*s, we work at O (mom 1 ). 
Any terms involving two-point, tree level vertices decorated by ^4*s vanish, at the 
desired order in momentum. Consequently, irrespective of the value of q, gauge 
invariance — in conjunction with what we have just done — ensures that the flow of 
our vertex vanishes at O (mom 1 ), without the need for any further tuning. Similarly, 
this result implies that gauge invariance can be used to demonstrate that the flow 
of a vertex decorated by four A l s and any number of C*s vanishes, without the need 
for further tuning. 

Our final task is to extend this analysis to include fermionic fields. To do this, 
we will treat -Bs and Ds separately, due to their differing mass dimensions. The 
point here is that neither S ^(k) nor S DD (k) vanishes at zero momentum. Thus, 
whereas gauge invariance can force some or all of the highest loop order seed action 
vertices with the maximum number of legs to be written as derivatives of lower 
point vertices when A l s are amongst the decorative fields, no such thing happens 
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here. Thus, both Bs and Ds behave essentially like the C*s, of the previous analysis. 

As a final point, we might worry about the vertex S ^ D (k) — which vanishes at 
zero momentum. However, supposing that it is the B(D) that is the internal field, 
this vertex will always be accompanied by a term in which there is a D(B) as an 
internal field. It is the seed action vertex at the other end of the corresponding 
dumbbell which is the one we tune. 

We have thus demonstrated that all dimensionless couplings, other than g and 
a, can be prevented from running by a suitable choice of the seed action. 
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4 Numerical Evaluation of fa 



4.1 The contributing Diagrams 



In this section, we give an expression for /?2 in terms of the A-derivative and a-terms 
which survive in the e — > limit. The complete set of diagrams contributing to fa 
(see [16]) can be derived using the techniques of [19], together with the subtraction 
techniques. Equivalently, it can be derived directly as a special case of the formula 
for arbitrary j3 n [16,20]. The expression for is: 



-WV(p) + 0(e) 



where (p) +2?2^u/(p) are gi ven i n figs- 21 and 22. Any diagrams contributing to 
02 are labelled T.#. 




T.4 T.5 T.6 T.7 



+ 




Figure 21: Diagrams contributing to T>2uv(p)- 
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4.2 The Universal Diagram 



Diagram T.8, when struck by — Ad\\ a , is the sole A-derivative term which yields 
simply a finite, universal contribution to /?2- 9 We reproduce this diagram, having 
chosen a particular momentum routing, in fig. 23. 




Figure 23: Reproduction of diagram T.8. 



The requirement that we take contributions which survive in the e — ► limit 
places useful constraints on the diagram. First, all fields must be in the ^-sector; 
given this, we are compelled to take O (mom ) from each of the vertices: additional 
powers of momentum would cause the diagram to vanish. For example, taking 0(p 2 ) 
from the vertices leaves us with, schematically 



0(p 2 



1 



l,k HI ~ k) 2 k* 



the integral is IR finite and is killed by Ad\\ a in the e — > limit. In turn, this forces 
both vertices to comprise a single supertrace: it is forbidden to have a single gauge 
field on a supertrace; if we take two supertraces, each with two gauge fields, then 
gauge invariance demands that there is no O(mom ) contribution to such a vertex. 

Now that we know that both vertices have only a single supertrace, all fields 
are forced to be in the A 1 sector. Temporarily ignoring attachment corrections, the 
group theory factor of the diagram must be either N 2 or unity. However, we can 
show that contributions of the latter type cancel. To see this, we can use the Ward 
identities to straightforwardly demonstrate [13, 16] 

(2) = —2(25^5^ — Sapd^ - 5^5^) 
Sv<y/3a(Q) = ~~ ^ (25 v p 5 ai — 5 1 p5 va — 5 vl 5p a ). 



9 Since the leading order contribution to this diagram is finite, it is not merely computable but 
actually universal: it is independent of the way in which we compute it (see [14] for a different way 
of evaluating essentially the same diagram). 
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Focusing on the component of diagram T.8 with a group theory factor of unity, 
the locations of the external fields are independent, since they are always guaranteed 
to be on the same supertrace. Summing over all independent locations of the external 
fields yields something proportional to 

S£j£ (0) + Sl 1 ^ (0) + (0) = 0. (40) 

Similarly, all attachment corrections can be ignored. If we suppose that one of 
the effective propagators attaches via a correction (see fig. 2) then the supertrace 
structure of the diagram is left invariant under independently placing the ends of 
this effective propagator in all independent locations. Hence the diagram vanishes 
courtesy of (40). Increasing the number of effective propagators which attach via a 
correction clearly does not change this result. 

Returning to the case of direct attachment, if the group theory goes as N 2 , then 
the locations of the external gauge fields are dependent, since it must be ensured 
that they are on the same supertrace. Up to insertions of -A* , we can use charge 
conjugation invariance to fix the order of the three internal fields so long as we 
multiply by two. Now, there are three identical pairs of locations that we can place 
the pair of fields A^ u . Including the diagram's overall factor of —1/6 we have: 



-N- 



x ^(0)5^(0) = -72N 2 V + 0(e). (41) 

To obtain the contribution to fa coming from diagram T.8, we must multiply the 
above factor by the number obtained from the loop integral. Since the integral yields 
a finite contribution, we simply Taylor expand the effective propagator A 11 (I — p) 
to 0(p 2 ). Remembering to evaluate the cutoff functions at zero momentum — which 
yields a factor of 1/2 for each of the effective propagators — we have: 



1 



4(L»V 

p.p 



k 2 {l - /c) 2 / 4 V' 1 I 



where we will define precisely what we mean by f, k in a moment. Looking at this 
expression, we might worry that the presence of Z 4 in the denominator means that 
the integral is actually IR divergent, even after differentiation with respect to Ad\\ a . 
However, due to the form of the 0(p 2 ) contributions, averaging over angles in the 
Z-integral will produce a factor of 

4 

in addition to the power of e coming from the A-derivative. This renders the contri- 
bution from diagram T.8 finite. 
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To evaluate the integral, we use the techniques of sec. 3. Specifically, we perform 
the I integral first, with unrestricted range of integration, and then perform the k- 
integral with the radial integral cutoff at A. After differentiation with respect to 
A5a| Q) the integral gives $%p 2 /32. Combining this with the factor coming from (41) 
yields: 



Before moving on, it is worth commenting further on the fact that all attachment 
corrections in diagram T.8 effectively vanish. When we finally come to evaluate the 
numerical value of fa, we will be dealing with diagrams for which all fields are in 
the A 1 sector. The highest point vertex that we will encounter is four-point: we 
have already seen how attachment corrections to such a vertex vanish. Three-point 
vertices are even easier to treat. Suppose that an effective propagator attaches 
via a correction to a three-point vertex, decorated exclusively by ^4s. We can sum 
over the two locations to which the effective propagator can attach, but these two 
contributions cancel, by charge conjugation. 

If nested gauge remainders are in the ^-sector, we know from [1, 16] that we can 
ignore attachment corrections. 

Thus, when we come to extract numerical contributions to @2, we will neglect 
attachment corrections. Similarly, for direct attachments, we need focus only on the 
cases where the group theory goes as A^ 2 . 

4.3 Finite, Non-Universal Diagrams 

Diagrams T.9 and T.10 both yield finite, non-universal contributions, when struck 
by —Ad\\ a . All contributions which survive in the D — > 4 limit must pick up a 
1/e pole (before differentiation with respect to A| Q ), which must come from the 
loop carrying the external momentum. Thus, we can take this loop to be in the 
A 1 sector and can Taylor expand the vertices which form this loop to lowest order 
in momenta. The combination of the bottom vertex and the effective propagators 
attached to it yield universal factors. However, it not helpful to convert these terms 
into algebra: keeping the diagrams intact will enable us to perform cancellations 
at the diagrammatic level. Indeed, it is only useful to explicitly Taylor expand the 
vertices (carrying p) which give a non-universal contribution, as shown in fig. 24. As 
usual, we have used charge conjugation to collect terms [1, 16]. 

The pattern of diagrams produced is highly suggestive: focusing in the top-most 
sub-diagrams, we can recognise terms that contribute to momentum derivatives of 
T>j . The missing terms are buried within the more complex, IR divergent diagrams 
contributing to fa, which we now examine. 




(42) 
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T.20 T.21 




T.22 T.23 




T.24 T.25 




Figure 24: Manipulation of diagrams T.9 and T.10 under Ad\\ a . This is valid up to 
0(e) corrections. 
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4.4 Subtractions 

Diagrams T.1-T.7, T.11-T.19 are all IR divergent, even after differentiation with 
respect to A\ a . Our strategy, as outlined in sec. 3.4, is to construct subtractions. 
Specifically, we construct subtractions for the latter set of terms, noting that this 
has already been done for diagram T.ll (see fig. 13). Our choice of subtractions for 
diagrams T.12-T.19 which, as noted already is not unique, is given in figure 29. 

Following sec. 3.4.2, we use the subtractions and additions to isolate the com- 
putable and non-computable parts of the parent diagrams. The latter terms are then 
manipulated, using eqn. (36). Processing any gauge remainders, further progress can 
be made by recognising that components of certain terms under the influence of NC 
vanish in the e — ► limit. As an example of this, consider the diagram shown on the 
left-hand side of fig. 25, which is obtained from the non-factorisable, non-computable 
component of diagram T.13. 




Figure 25: Manipulation of a term under the influence of NC. 



To understand how to manipulate the diagram on the left-hand side, it suffices 
to look at the component Taylor expandable in p. We begin by temporarily ignoring 
the NC and looking at the most IR divergent contribution. In this case, the integrand 
goes like: 

l\l-k) 2 k A {P ' 

The NC forces additional powers of l,k to appear in the numerator. The 0(1 ■ k) 
now does not survive differentiation with respect to A| Q in the e — > limit and 
the 0(l 2 ) term survives if we take additional powers of I, only. Consequently, the I 
integral becomes Taylor expandable in k, yielding the right-hand side of fig. 25. This 
manipulation has yielded a term which can be processed using (36). 
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Iterating the diagrammatic procedure yields: 
1 d _ , x 9N 2 



-4&CV(p)+0(e) = -271^1^(P)-^PV+[^(P) C + P 2>(P)J > ( 43 ) 
where 'D^aviP) an d ^2u,u(p) are given, respectively, in figs. 26 and 27. 



It is thus apparent that (a-terms aside) the expression for 02 is computable, up 
to diagrams T.42-T.50. If we are to obtain a universal (in the a — > limit) 
it must therefore be the case that these latter diagrams are computable also. To 
demonstrate this, we need not construct further subtractions but rather need only 
utilise the transversality of ~~ |_JJ — when the external legs are in the ^-sector — as 

they are for all contributions to I'lui/W which survive in the e — ► limit. This 
transversality can be exploited to extract the momentum dependence and Lorentz 
structure of 



1 



1 



1 



Ql-+Q2 + (H1- + H2 



2c 



2U af3 (q) + ... 



The ellipsis denotes terms higher order in momentum and / or e, and the factor of 
two is extracted for convenience. Notice that this transversality immediately reduces 
the apparent severity of the IR divergences of diagrams T.42-T.44 (c/. the comments 
above (31)). 

For all non-factorisable diagrams, it must be the case that ~~ |_JJ — attaches to 
something computable, else the diagram as a whole vanishes in the e — ► limit. 
This is easiest to see by working in D = 4, where the potentially most IR divergent 
contributions to such diagrams are obtained by taking the part of lJJ — which 

goes as 0(/ 2 )lnL Focusing, for example, on the part of such a diagram Taylor 
expandable in p, we will only find an IR divergence after integrating over I if the 
/-integrand picks up a factor of 1/Z 6 ; such contributions are computable. 

With this in mind, and recalling that Q\ and the Hi are computable, the non- 
universal contributions to diagrams T.42-T.50, which survive in the e — ► limit are 
as follows: 

1. the components of diagrams T.46-T.50 in which we take the non-computable 
part of the bottom most sub-diagram; 

2. the components of diagrams T.42-T.49 in which we take the G2 contribution 
and take only the computable parts of all remaining diagrammatic 



1 



to 



elements. 
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T.26 



T.27 



T.28 



T.29 
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T.42 



T.43 



T.44 
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In both cases, the sum of all contributions can be shown to vanish, diagram- 
matically! In the first case, we note that the non-computable components of the 
bottom-most sub-diagrams of T.46-T.47 are Taylor expandable in p (the same can- 
not, in general, be said of the computable components). Performing the usual dia- 
grammatic manipulations, it is straightforward to show that the elements of item 1 
cancel, amongst themselves. 

To show that the diagrams of item 2 cancel, we exploit the transversality of 



1J noting that, if \ 1 | is contracted into an effective propagator, we can 



use the effective propagator relation, up to terms which we do not care about: 

2n a p(q)A n (q) = l + 0(q 2 ), 

where, if we identify q with a loop momentum, the 0(q 2 ) terms do not contribute 
to @2 in the e — > limit whereas, if we identify q with the external momentum, the 
0(q 2 ) terms yield diagrams which vanish at 0(p 2 ). 

Applying the effective propagator relation in diagrams T.46-T.49, the resulting 
gauge remainder can be discarded: the sub-diagram it strikes is transparent to it, by 
Lorentz invariance, and so it effectively hits the 0{p 2 ) stub, killing it. Applying the 
effective propagator relation in diagram T.43, the Q2 part of this diagram cancels 
those of T.48 and T.49, up to corrections which vanish in the e — > limit. Finally, 
we turn to the Q2 part of diagram T.42. Applying the effective propagator relation, 
the Kronecker-<5 contribution cancels the Q2 parts of diagrams T.44 and T.45 at 
leading order in e. Processing the gauge remainders removes the Q2 contributions 
from diagrams T.46 and T.47 in the e — > limit. 

We have thus demonstrated that, up to the a-terms, fa can be reduced to a 
sum of computable contributions. The IR divergent terms and terms which go as 
p~ 2e (see sec. 3.3) cancel between diagrams. The following finite contributions to 
— 4:02^ nu(j>) are left over, and can be evaluated using the techniques of sec. 3.4: 



diagram(s) 


N 2 /(Att) 4 


T.8 




T.26-T.41 


-733/9D^(p)+9p 2 V 


T.42-T.49 


1141/9CV(p) 



Putting everything together, we have: 



34 N 2 Id 
/^□^(p) + 0(e) = -Y7^j4 D pW + 8 J1 d^ Vl ^ p) - (44) 
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4.5 The a-terms 
4.5.1 The problem 

To deal with the a-terms, we must understand the a-dependence of Vi^ip). To 
make this procedure as transparent as possible, we will use the explicit forms for 
the two-point, tree level vertices, effective propagators and gauge remainders given 
in appendix A. We note, though, that it should be possible to repeat this analy- 
sis, using only the general properties which we know these functions must satisfy 
as a consequence of Lorentz invariance, gauge invariance, dimensions and the UV 
finiteness of the theory. 

In the current picture, we have an explicit algebraic realisation for four of the 
members of T>i fiu (p): the final member of the standard set and all of the members 
of the little set. The remaining diagrams all contain either three-point or four-point 
vertices. It is not our aim to choose specific forms for these vertices through a specific 
choice of seed action and boundary conditions for the flow; rather, our aim is to show 
that only implicit choices are necessary for our purposes. 

Assessing the a-dependence of T>i^ u (p) is complicated by the fact that a loop 
integral must be performed. This integral must be done before we take the a — > 
limit, as the two procedures do not commute. However, whilst we cannot set a — ► 
too soon, we can work at small a, and will do so henceforth. In this limit, it follows 
from eqn. (15) that 71 ~ a. Consequently, the a-terms will vanish, as required, only 
if the behaviour of dT>i^ u (p)/da is better than 1/a. 

From our work in sec. 3, we know how to parameterise the a-dependence of 
T>i^ u (p). For our purposes, this is most easily done in D = 4: 



The non-universal function H(a) is independent of A. We now choose to recast 
this equation. When constructing the (dimensionless) argument of the logarithm, 
we divide the IR scale by the only other scale available, A. 

Let us examine this in the context of actually performing the loop integrals to 
obtain (45). The appearance of the IR scale has been discussed, in depth, in sec. 3. 
The scale A has a natural interpretation as the scale at which the loop integrals 
are effectively cutoff in the UV. However, we should not preclude the possibility 
that the loop integrals are actually cutoff at some scale /i(a)A, where h{a) is a 
dimensionless function, independent of A. Of course, this has no effect on the value 
of (3\ obtained by differentiating T>\ tiu {p) with respect to A| Q . With this in mind, we 




(45) 
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rewrite eqn. (45) as follows: 

□/«/(?)• (46) 

This recasting now allows us to break the problem of the a-terms into two parts. 
On the one hand, we have potential a-dependence coming from any non-trivial in- 
dependence of the effective cutoff scale, parameterised by h(a). On the other hand, 
we have a-dependence coming from the region of the loop integral with support, 
parameterised by H(a). We deal with these cases in turn. 



4/3i In 



(IR scale)' 
A 2 h{a) 



+ H(a) 



4.5.2 Behaviour of h(a) 

The treatment of this problem is slightly easier than one might expect. The crucial 
point is that the logarithm term in eqn. (46) comes only from (UV regularised) 
terms with non-trivial IR behaviour. There are five diagrams with non-trivial IR 
behaviour: the final two elements of the standard set and the elements of the little 
set. Our strategy is to examine these diagrams and, through a choice of cutoff 
functions, ensure that the momentum integrals are cutoff at A; equivalently that 
h(a) is independent of a. 

Sufficient UV regularisation can be provided by cutoff function regularisation 
alone — e.g. for the little set — or entirely by the regularising sector — e.g. for the final 
element of the standard set. In the latter case, we are interested simply in the scale 
at which the .B-sector diagrams regularise the ^-sector diagrams. In the former 
case, we are interested not only in this but also the scale at which the A 1 sector 
diagrams die off on their own. 

For all that follows, we assume that the momentum of the cutoff functions 
c(k 2 /A 2 ) and c(k 2 /A 2 ) crosses over from large to small for x = k 2 / A 2 ~ 0(1)- 
This amounts to an implicit choice of the non- universal details of the set-up. In this 
section, we will demonstrate that, given this choice, we can consistently arrange for 
all momentum integrals to be cutoff at x ~ 0(1). 

We start by looking at the final element of the standard set which has the alge- 
braic form 

i N f f ( k ~ p ^ k t> f k ~p( k ~ Ek f kk ^ \ (A7\ 
Jk \{k-p) 2 k 2 A 2 A 2 J' { ' 

where 

= (1 + a)cx , 4g , 
(1 + a)xc x + 4ac x 

At large x, where fk = 1/x [1, 13], we recover unbroken SU(N\N), as we must, for the 
theory to be regularised. Given our assumption that c x crosses over at x ~ 0(1), we 
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need only ensure that the denominator of fk crosses over at x ~ 0(1). The crossover 
of the denominator occurs at 

xc x ~ 4ac x . 

Now, since we are working at small a (and xc x ^> c x for large x), the crossover must 
happen for small values of x. Taylor expanding, we therefore find that the crossover 
occurs at 

xcq ~ aco, 

(in the limit of small a), cq is fixed to be unity by the renormalisation condition 
for A 1 , eqn. (5); however, there is no such constraint on cq. In turn, this implies 
that the momentum integral for the third element of the standard set is cutoff at 
x ~ O(a/co). Note that if we set Co ~ 1, then we would indeed find the problem 
that dT>i /lu (p)/da ~ 1/a. Demanding that the cutoff scale occurs at x ~ 0(1) then 
forces us to choose 

c ~ 0(a) (49) 

(which is perfectly compatible with c x crossing over at x ~ 0(1)). 

Let us now turn to the remaining diagrams with non-trivial IR behaviour. The 
treatment of these is somewhat different from what we have just done, as we have 
regularisation provided not only by the -B-sector, but also by cutoff function regu- 
larisation. The crossover scale in the -B-sector follows trivially from the observation 
that 

y ' (l + a)A 2 

Hence, we immediately know that, given the choice cq ~ 0(a), the crossover occurs 
at x ~ 0(1), in this sector. 

However, we now need to show that the scale at which the cutoff regularisa- 
tion kicks in in the A 1 -sector also occurs at this scale. If a diagram is sufficiently 
regularised by cutoff regularisation alone, then the B-sector becomes effectively re- 
dundant. If the B-sector is required, in addition to cutoff regularisation, then there 
can be two scales in the problem: the first is where the momentum in the ^-sector 
can be considered large and the second is where the momentum in the -B-sector can 
be considered large. The A 1 and B-sectors cancel each other at the highest of these 
scales. 

Turning now to the ^-sector, 

1 ' k^(a + l) + c x (a-iy 

which goes as ac x /x for large x and as l/2x for small x. The crossover occurs, for 
small a, at 

a(c x + 1) ~ 1 - c x . 
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(Note that the left-hand side dominates at sufficiently small x, whereas the right- 
hand side dominates at large x.) Again, due to the smallness of a, we can Taylor 
expand in x to find the crossover point, which occurs at 



X |Cn 



o 



0(1). 



This implies that, for the effective cutoff to be at x ~ 0(1), we must choose \c' \ ~ 
0(a). 

This actually completes the analysis necessary to show that h{a) can always be 
arranged to be independent of a. However, for the purposes of the next section, it 
is useful to show that we can, in fact, ensure that all momentum integrals are cutoff 
at x ~ 0(1)- The reason that this is useful is because we expect a one dimensional 
integral with an integrand of 0(1) but with support only over a range Ax to go like 
Ax. Hence, it is desirable for this range to be 0(1) as opposed to, e.g. , (a" 1 ). 

In the A 2 sector, from 



22/- ac * 



A 22 (k) 



k 2 a + 1 + c x (l — a) 



we see that the large momentum behaviour is ac x /x, whereas the small momentum 
behaviour is a/2x. Since the a-dependence is the same for both, the crossover scale 
is clearly set by c, which is assumed to crossover at x ~ 0(1). 
In the D sector, given that 



DD(i„\ _ &xfx 



A^(fc) 



A 4 



our choice that cq ~ 0(a) ensures that the crossover occurs at x ~ 0(1), assuming 
that c x crosses over at x ~ 0(1). 
In the C*-sector, 

iV X -j- LAL-x 

from which it is clear that the crossover scale is controlled by c and A; we use this 
freedom to ensure that the crossover occurs at x ~ 0(1). 

We have thus demonstrated that, by suitable choices of the behaviours of our 
cutoff functions, we can guarantee that all momentum integrals are cutoff at the scale 
x ~ 0(1) (working at small a); one consequence of this is that h(a) is independent 
of a and so does not generate a contribution to /?2 in the a — > limit. 

We conclude this section with an interesting comment on universality. It is clear 
from our analysis so far that the freedom to choose the non-universal parts of our 
cutoff functions enables us to choose h(a). Returning to eqn. (46), it thus looks like 
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we could generate a universal contribution to P2 by choosing e.g. h(a) = a m , for 
some m ^ 0. However, the universal appearance of this contribution is accidental, 
as can be appreciated from the fact that it arises from a particular choice of a non- 
universal function. Indeed, the universal 02 will only be obtained by arranging things 
so that all contributions from the running of a can be removed in the a — > limit. 

4.5.3 Behaviour of H(a) 

To start our analysis of H(a), we begin by returning to the third element of the 
standard set. We know that the integrand of eqn. (47) effectively has support only 
over the region < x < 0(1). Moreover, any non-trivial a-dependence of H{a) must 
come from the -B-sector, as ^4-sector (processed) gauge remainders are independent 
of a. 

From our algebraic choice for / (see equation (48)), the most obvious possible 
source of problematic a-behaviour comes when x is small. However, this is amelio- 
rated by our previous choice of cq. To complete our analysis of this diagram, we must 
now perform the loop integral. However, since the effective cutoff of this integral is 
O(l) as opposed to, say, (a^ 1 ), we do not expect the loop integral to generate any 
bad a-dependence (i.e. dependence which diverges as a — > 0). 

Looking now at the little set, the situation is similar: in the -B-sector, our choice 
of Co cures any bad a-dependence; in the ^-sector there is not even a potential 
problem. In both cases, the effective cutoff for the loop integral is 0(1). This 
exhausts the analysis of the diagrams for which we have an explicit algebraic form 
and so now we turn to the diagrams possessing three and four-point vertices. 

First, we will look at the diagrams with an A^A^C 1 vertex. Neither this vertex 
nor the effective propagator to which it attaches carries the loop momentum of the 
diagram and so we analyse them separately. The effective propagator — which carries 
zero momentum — goes as 

1 

2AA 4 ' 

We see that the a-dependence of this can always be controlled by a suitable choice 
of A. (It can always be ensured that this choice of A is compatible with the choice 
which guarantees that the crossover for the C*-sector effective propagator occurs at 
x ~ 0(1).) 

What about the A^A^C 1 vertex? We have encountered this already in section 3.5: 
at 0{p 2 ) it is a dimensionless coupling and so we tune its flow to zero. This means 
that the flow equations tell us nothing about its a-dependence: this dependence is 
a boundary condition. The solution is simply to choose the boundary condition to 
have sufficiently good a-dependence, which is something we are always at liberty to 
do. 
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Now let us examine the remaining part of these diagrams. Attached to the other 
end of the C l -sector effective propagator is either a hook or a three-point vertex, 
decorated by a simple loop. Since the hook simply goes as 



where g x = (1 — xf)/2, it is clear that the loop integral does not produce any 
troublesome a-dependence. 

The case where the top part of the diagram constitutes a three-point vertex is 
almost as easy to treat. Let us consider the flow of this vertex. If we take the 
dimensionless part, then we know that the flow has already been tuned to zero; 
in this case, we choose the boundary condition, appropriately. Taking the flow of 
the dimensionful part of the vertex we once again tune the seed action. This time, 
though, we do so to ensure sufficiently good a-dependence. Note that we need not 
worry about constraints coming from gauge invariance. If the top vertex contains 
two fields in the ^4-sector (it cannot contain a single one) then gauge invariance 
simply tells us that the vertex is transverse; it is not related to lower point vertices. 
Performing the loop integral does not generate any bad a-dependence. 

The penultimate diagram to deal with is the second element of the standard set, 
which comprises two three-point vertices, joined together by two effective propaga- 
tors. Once again, our aim is to choose the seed action such that the a-dependence 
of the three-point vertices is sufficiently good. This time, however, we must worry 
about gauge invariance. 

The first effect of gauge invariance is to relate the longitudinal part of each of 
these vertices to two-point, tree level vertices. Referring to our list of two-point, tree 
level vertices (eqns. (50)-(55)), it is clear that our three-point vertices have buried 
in them, necessarily, components which go as O (a" 1 ). 

To examine the transverse components of these vertices, we must use the flow 
equations. Fig. 28 shows the flow of a three-point vertex comprising one external 
field and two identical internal fields (i.e. this vertex should be viewed as part of a 
whole diagram). 



We expect the critical case to occur when all three fields are in the ^-sector, 
since gauge invariance will then force us to take O (mom 2 ) from each two-point, tree 
level vertex. In particular, this means that we will be unable to tune the three-point 
seed action vertices if we take the O (mom 3 ) part of the vertex whose flow we are 
computing (the O (mom 1 ) part is, of course, universal and independent of a). 
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Figure 28: Flow of a three-point, tree level vertex viewed as part of a whole diagram. 
Given our choice that c' ~ 0(a), it is straightforward to show that 



Thus it is clear that we can always tune the seed action to ensure that worst 
behaviour of the three-point, tree level vertices is the O (a -1 ) dependence forced by 
gauge invariance. This leading a-dependence now cancels between the vertices and 
effective propagators of the second element of the standard set. The loop integral 
does not generate any bad a-dependence 

The treatment of the first element of the standard set follows, similarly: by 
tuning the seed action (and choosing suitable boundary conditions for the flow) we 
can ensure that the worst a-dependence is that forced on us by gauge invariance. 
This dependence is then cancelled by the effective propagator. The loop integral 
does not generate any bad a-dependence. 

We have thus demonstrated that, by suitable choice of seed action, boundary 
conditions for our flow and non- universal behaviour of the cutoff functions, we can 
ensure that H(a) ~ O (a ). It therefore follows that we can always arrange for the 
a-terms to vanish in the a — > limit. Finally, then, we recover the standard value 



5 Conclusions 

We have performed the first manifestly gauge invariant, continuum computation 
of the SU{N) Yang-Mills two-loop function; no gauge fixing or ghosts being in- 
troduced at any stage. The framework employed is the ERG of [1, 16] which was 
constructed to allow convenient renormalisation beyond one loop. This formalism is 




O (a ) O (mom 3 ) 
A2 



of ft: 
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furnished with a set of powerful diagrammatic techniques, reviewed in sec. 2. These 
techniques form the basis of the methodology used to reduce both (3\ and /?2 to a 
set of terms from which the universal coefficient can be readily extracted. Let us 
recapitulate the basic procedure at one loop. 

From the ERG equation, a diagrammatic expression is generated for the flow of 
the classical, two-point, A 1 vertex. Taking this vertex to carry momentum p, we spe- 
cialise to one loop and use the renormalisation condition. This yields a diagrammatic 
expression for (3\ which, since it contains instances of the seed action and details of 
the covariantisation, is not yet manifestly universal. To proceed, we convert terms 
comprising exclusively Wilsonian effective action vertices and undecorated kernels 
into A-derivative terms and corrections. A subset of the correction terms can be 
simplified using the effective propagator relation, cancelling non-universal terms, up 
to gauge remainders. In turn, these remainders can be processed diagrammatically. 
Iterating the procedure, the expression for (5\ can be reduced to A-derivative terms 
and terms possessing an 0(p 2 ) stub. At one-loop, these latter terms are treated by 
Taylor expanding the sub-diagram attached to the stub to zeroth order in p, thereby 
allowing the 0(p 2 ) terms to be reduced to A-derivative terms also. 

At two loops, the procedure is much the same [1, 19], though it is complicated 
by the fact that direct Taylor expansion of all 0{p 2 ) terms is no longer possible. In 
certain diagrams, p effectively acts as an IR regulator, and so Taylor expansion in p 
generates spurious IR divergences. Isolation of the non- Taylor expandable compo- 
nents is most easily achieved by use of subtraction techniques, which are introduced 
in sec. 3.4. This then allows to be reduced to a set of A-derivative and a-terms [16]. 

The strategy for evaluating A-derivative terms is to start by interchanging the 
order of differentiation with respect to A| Q and loop integration. This makes it clear 
that, so long as there are no dimensionless running couplings hidden in the integrand, 
only those terms which, before differentiation with respect to A| a , possess some IR 
scale will survive. In sec. 3.5 it is proven that all such running couplings can be 
removed from the set-up. 

At one loop, the extraction of the numerical coefficient is straightforward. Indi- 
vidual terms either vanish or contribute a finite, universal number, the sum of which 
combine to yield the standard answer. At two loops, extraction of a numerical coef- 
ficient from the A-derivative terms is much harder, since individual terms can now 
develop IR divergences which survive even after differentiation with respect to A| Q . 
Whilst it must be the case that the sum of diagrams is IR finite, the sub-leading 
terms are no longer manifestly universal. 

Though it is not possible to unambiguously define the universal (finite) com- 
ponent of a diagram which also makes a finite, non-universal contribution to fo, 
progress can be made. Utilising the subtraction techniques, a prescription — which 
is by no means unique — can be defined which separates a diagram into computable 



64 



and non-computable parts. By non-computable, it is meant that the corresponding 
coefficient cannot be calculated without specifying the appropriate parts of the seed 
action and details of the covariantisation. The computable part can be evaluated 
using only the renormalisation condition for A 1 , but depends on the precise prescrip- 
tion used to define what is meant by computable. Separating terms in this way, it 
can be shown, diagrammatically, that all non-computable contributions to cancel, 
amongst themselves. The sum of the computable contributions yields the standard, 
universal answer. 

The final task is to treat the a-terms. These terms are expected to be present as 
a consequence of an unphysical coupling associated with one of the regulator fields. 
Since fa is not a physical quantity, it is anticipated that agreement between the value 
obtained in our scheme and the standard value is found only in the limit that a is 
tuned to zero [1, 16]. In sec. 4.5 it is demonstrated that, by sufficiently constraining 
the seed action and boundary conditions of the flow, the a-terms can indeed be 
shown to vanish in this limit, so long as we impose constraints on the non-universal 
behaviour of the cutoff functions. 

Thus, we have successfully computed /?2 without fixing the gauge, finding agree- 
ment with the standard, universal answer in the appropriate limit. This demon- 
strates the consistency of the formalism beyond reasonable doubt, thereby opening 
up the prospect of an analytical, manifestly gauge invariant computational scheme 
for SU(N) Yang-Mills theory. In the future, we aim to apply the formalism both 
perturbatively and non-perturbatively and also plan to incorporate quarks. 
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A Ingredients of the Weak Coupling Flow Equations 

A list of our choices for all single supertrace seed action, two-point, tree level vertices 
follows. Multi-supertrace terms are either related to those listed, by no- .A symmetry, 
or can be set to zero [16]. The undetermined parameter A can depend on a. 




ac. 



(50) 



(51) 
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s, b J(p) 



a + 1 



>=W(p) + — 
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S Q D ?(p) = (53) 

Cp 

S DD (P) = ^f, (54) 
S cic \p) = ^ + 2AA 4 . (55) 

Cp 

These choices correspond to [16] 

f- = nJV"^ ■ < 56 > 

(1 + ajxcp + Aacp 



2ac p 



(1 + a)xc p + 4ac p ' 
allowing us to write the effective propagators in the following forms: 



(57) 



A n (p) = \ rr, (58) 

x ' p 2 a + 1 + c p {a - 1) y J 

A 22 (p) = 1 (59) 

w a + 1 + Cp(l - a) y ' 

& BS (P) = ^p9 p , (60) 

A D5 (P) = (61) 

A 4 x + 2Ac p 



B Subtractions for diagrams T.12— T.19 
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